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1 Hadron Physics 

Detailed investigations of the structure of hadrons are essential for under- 
standing how matter is constructed from the quarks and gluons of Quantum 
chromodynamics (QCD), and amongst the questions posed to modern hadron 
physics, three stand out. What is the rigorous, quantitative mechanism re- 
sponsible for confinement? What is the connection between confinement and 
dynamical chiral symmetry breaking? And are these phenomena together suf- 
ficient to explain the origin of more than 98% of the mass of the observable 
universe? Such questions may only be answered using the full machinery of 
nonperturbative relativistic quantum field theory. 

This contribution provides a perspective on progress toward answering 
these key questions. In so doing it will provide an overview of the contempo- 
rary application of Dyson-Schwingcr equations in Hadron Physics, additional 
information on which may be found in Refs. [1, 2, 3, 4, 5, 6]. The presentation 
assumes that the reader is familiar with the concepts and notation of relativis- 
tic quantum mechanics, with the functional integral formulation of quantum 
field theory and with regularisation and renormalisation in its perturbative 
formulation. For these topics, in order of appearance, Refs. [7, 8, 9, 10] are 
useful. In addition, Chaps. 1 and 2 of Ref. [5] review the bulk of the necessary 
concepts. 

Hadron physics is a key part of the international effort in basic science. For 
example, in the USA we currently have the Thomas Jefferson National Accel- 
erator Facility (JLab) and the Relativistic Heavy Ion Collider (RHIC) while 
in Europe hadron physics is studied at the Frascati National Laboratory and 
is an important part of a forthcoming pan-European initiative; namely, the 
Facility for Antiproton and Ion Research (FAIR) at GSI-Darmstadt. Progress 
in this field is gauged via the successful completion of precision measurements 
of fundamental properties of hadrons; e.g., the pion, proton and neutron, and 
simple nuclei, for comparison with theoretical calculations to provide a quan- 
titative understanding of their quark substructure. 
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The proton and neutron (collectively termed nucleons) are fermions. They 
are characterised by two static properties: an electric charge and a magnetic 
moment. In Dirac's theory of pointlike relativistic fermions the magnetic mo- 
ment is 

"° = m- (L1) 

The proton's magnetic moment was discovered in 1933 by Otto Stern, who 
was awarded the Nobel Prize in 1943 for this work [11]. He found, however, 
that 

M P = (l + L79)/i B . (1.2) 

This was the first indication that the proton is not a point particle. Of course, 
the proton is now explained as a bound state of quarks and gluons, a com- 
position which is seeded and determined by the properties of three valence 
quarks. These quarks and gluons are the elementary quanta of QCD. 

Important hadron physics experiments measure hadron and nuclear form 
factors via electron scattering. The electron current is known from QED: 

j fl (P',P) = ieu e (P')A li (Q,P)u e (P), Q = P'-P (1.3) 
= ieu e (P') lfl (-l)u e (P), (1.4) 

where u e , u e are Dirac spinors for a real (on-mass-shell) electron. Equation 
(1.4) is the Born approximation result, in which the dressed-electron-photon 
vertex is just (— 7 M ) and where the negative sign merely indicates the sign of 
the electron charge. For the nucleons, Eq. (1.3) is written 

jJf(P',P)=ieu N (P')A^Q,P)u N (P) (1.5) 

= ieu N (P>) (i^(Q 2 ) + -La^Q„F 2 N (Q 2 )^ u N (P) (1.6) 

where un, um are on-shell nucleon spinors, and ' (Q 2 ) is the Dirac form 
factor and F 2 N (Q 2 ) is the Pauli form factor. The so-called Sachs form factors 
are defined via 

G E (Q 2 ) = F X (Q 2 ) - ^F 2 (Q 2 ) , G M (Q 2 ) = F 1 (Q 2 ) + F 2 (Q 2 ) . (1.7) 

In the Brcit frame in the nonrclativistic limit, the three-dimensional Fourier 
transform of Ge(Q 2 ) provides the electric-chargc-density distribution within 
the nucleon, while that of Gm{Q 2 ) gives the magnetic-current-density dis- 
tribution. This explains their names: G^(Q 2 ) is the nucleon's electric form 
factor and G^j(Q 2 ) is the magnetic form factor. It is apparent via a compar- 
ison between Eqs. (1.4) and (1.6) that F 2 = for a point particle, in which 
case Ge — Gm- This means, of course, that if the neutron is a point particle 
then it possesses neither a charge nor a magnetic moment; i.e., they are both 
zero. 
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There arc six quarks in the Standard Model of particle physics: u (up), 
d (down), s (strange), c (charm), b (bottom) and t (top). The first three are 
most important in hadron physics. A central goal of nuclear physics is to 
understand the structure and properties of hadrons, and ultimately atomic 
nuclei, in terms of the quarks and gluons of QCD. So, why don't we just go 
ahead and do it? One of the answers is confinement: no quark or gluon has 
ever been seen in isolation. Another is dynamical chiral symmetry breaking 
(DCSB); e.g., the masses of the u, d quarks in perturbative QCD provide no 
explanation for ~ 98% of the proton's mass. One therefore has to ask, with 
quarks and gluons are we dealing with the right degrees of freedom? 

The search for patterns in the hadron spectrum adds emphasis to this ques- 
tion. Let's consider the proton, for example. It has a mass of approximately 
1 GeV. Suppose it to be composed of three (two u and one d) constituent- 
quarks, as in the eightfold way classification of hadrons into groups on the basis 
of their symmetry properties [12]. A first guess would place the mass of these 
constituents at ~ 350 McV. In the same approach, the pseudoscalar 7r-mcson 
is composed of a constituent-quark and a constitucnt-antiquark. It should 
therefore have a mass of ~ 700 McV. However, its true mass is ~ 140 McV! 
On the other hand, the mass of the vector p-meson is correctly estimated in 
this way: m p = 770 MeV. Such mismatches are repeated in the spectrum. 

Furthermore, modern high-luminosity experimental facilities that employ 
large momentum transfer reactions are providing remarkable and intriguing 
new information on nucleon structure [13, 14]. For an example one need only 
look so far as the discrepancy between the ratio of electromagnetic proton 
form factors, fi p G p E (Q 2 )/G p M (Q 2 ), extracted via the Rosenbluth separation 
method [15] and that inferred from polarisation transfer experiments [16, 17, 
18, 19, 20], Fig. 1. This discrepancy is marked for Q 2 > 2 GeV 2 and grows 
with increasing Q 2 . 3 Before the JLab data were analysed it was assumed that 

Cb( ( 3 2 )/ ( ^m( ( 3 2 ) = 1 based on the seemingly sensible argument that the 
distribution of quark charge and the distribution of quark current should be 
the same. However, if the discrepancy is truly resolved in favour of the JLab 
results, then we must dramatically rethink that picture. 

An immediate question to ask is where does the modern theory of hadron 
structure stand on this issue? A theoretical understanding might begin with a 
calculation of the proton's Poincare covariant wave function. (Remember, the 
discrepancy is marked at larger momentum transfer; viz., on the relativistic 
domain.) One might think that is not a problem. After all, the wave func- 
tions of few electron atoms can be calculated very reliably. However, there 



3 It is currently believed that the discrepancy between the two classes of exper- 
iment may be removed by the inclusion of two-photon exchange contributions 
in the analysis of the ep scattering cross-section; i.e., by improving on the Born 
approximation. Estimates indicate that such effects materially reduce the magni- 
tude of the ratio inferred from Rosenbluth separation and slightly increase that 
inferred from the polarisation transfer measurements. See, e.g., Ref. [21] 
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Fig. 1. Proton form factor ratio: fi p G P E (Q 2 ) / G p h[ (Q 2 ). The data are obtained 
via Rosenbhith separation, circles - Ref. [15]; and polarisation transfer, squares - 
Ref. [16]; diamonds - Ref. [18]. (Adapted from Ref. [5].) 



are some big differences. One is that the "potential" between light-quarks is 
completely unknown throughout ~ 98% of the proton's volume because in 
QCD perturbation theory is unreliable at length-scales > 0.2 fm. Moreover, as 
we shall see, a reliable description of the proton's wave function will require 
an accurate treatment of virtual particle effects, which are a core feature of 
relativistic quantum field theory. In fact, a computation of the proton's wave 
function requires the ab initio solution of a fully-hedged relativistic quan- 
tum held theory. That is yet far beyond the capacity of modern physics and 
mathematics. 



1.1 Aspects of QCD 

The theory one must explore is quantum chromodynamics. It is a local, renor- 
malisable, non-Abelian gauge theory, in which each flavour of quark comes in 
three colours and there are eight gauge bosons, called gluons. It is a pecu- 
liar feature of non-Abelian gauge theories that the gauge bosons each carry 
the gauge charge, colour in this case, and hence self-interact. This is the key 
difference between QED, an Abelian gauge theory wherein the photons are 
neutral, and QCD. The gluon self-interaction is primarily responsible for the 
marked difference between the running coupling in QCD and that in QED; 
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namely, that the running coupling in QCD decreases relatively rapidly with 
increasing momentum transfer - the theory is asymptotically free, whereas 
the QED coupling increases very slowly with growing momentum transfer. 
(See, e.g., Sect. 1.2 in Rcf. [22].) 

The QCD action is expressed through a local Lagrangian density; viz., 

S[A* q, q) = [ d A x \q{x) [ 7/1 £> M + M] q(x) 



+\f^ v (x)F^(x) +±d„A«(x)d 1/ AZ(x)}, (1.8) 

in which the first term of the second line is the chromomagnctic field strength 
tensor 

f; u (x) = dMx) - duA%{x) + gf abc Al(x)Al(x), (1.9) 

where {f abc : a,b,c — 1,...,8} are the structure constants of SU(3), and 
the second term is a covariant gauge fixing term, with £ the gauge fixing 
parameter: £ = 1 defines Feynman gauge and the limit £ — > implements 
Landau gauge in the generating functional. The first term in Eq. (1.8) involves 
the covariant derivative: 

I>,. = M -i0yA-(*), (1.10) 

where {|A a ;a = 1, ... ,8} are the generators of SU(3) for the fundamental 
representation; and, in addition, the current-quark mass matrix 



M 



fm u ...\ 
m d .. 
m s . . 

V ; j 



(l.n) 



wherein only the light-quark elements are indicated. 

Understanding hadron obscrvablcs means knowing all that the quantum 
field theory based on Eq. (1.8) predicts. As has been emphasised, perturba- 
tion theory is inadequate to that task because confinement, DCSB, and the 
formation and structure of bound states are all essentially nonperturbative 
phenomena. An overview of a nonperturbative approach to exploring strong 
QCD in the continuum will subsequently be provided. 4 

In Table 1 we return to the hadron spectrum and focus on some of its 
features. In a constituent-quark model the J PC = 1 p-meson is obtained 
from the + 7r-meson by a spin flip, yielding a vector meson state in which the 



4 Almost all nonperturbative studies in relativistic quantum field theory employ a 
Euclidean Metric. A Euclidean metric was used in writing Eq. (1.8). Appendix A 
provides some background and describes the Euclidean conventions employed 
herein. 
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Table 1. Comparison between various hadron mass ratios, involving the it, p, a and 
ai mesons, the nucleon - N, and the first radial excitations of: the tt = 7Ti; and the 
p = p\. The individual masses can be found in Ref. [23]. 

9 9 

• — £ = 30 • — 2i = 2.1 ? Hypcrfinc Splitting 

• — ^- =86 • — = 3.5 ? Excitation Energy 

• — — « 7 • — — = — ~ — ? Quark Counting 
m n m„ 4 2 



spins of the constituent-quark and -antiquark are aligned. The same procedure 
would yield the 1 ++ ai-meson from the ++ cr-mcson. Thus the difference 
between the masses of the p and 7r, and the a\ and a would appear to owe to 
an hyperfine interaction. As the Table asks in row one, why is this interaction 
so much greater in the it channel? Another question is raised in row two: why 
is the radial excitation energy in the pseudoscalar channel so much greater 
than that in the vector channel? And row three asks the question that was 
posed on page 3: why doesn't constituent-quark counting work for the pion's 
mass? 

Additional questions can be posed. The range of an interaction is inversely 
proportional to the mass of the boson that mediates the force. The nucleon- 
nucleon interaction has a long-range component generated by the it. The fact 
that the pion is so much lighter than all other hadrons composed of u- and 
d-quarks is crucially important in nuclear physics. If this were not the case; 
viz., were the pion roughly as massive as all like-constituted hadrons, then 
the domain of stable nuclei would be much reduced. In such a universe the 
Coulomb force would prevent the formation of elements like Fc, and planets 
such as ours and we, ourselves, would not exist. 

1.2 Emergent Phenomena 

A true understanding of the visible universe thus requires that we learn just 
what it is about QCD which enables the formation of an unnaturally light 
pseudoscalar meson from two rather massive constituents. The correct under- 
standing of hadron observables must explain why the pion is light but the 
p-meson and the nucleon are heavy. The keys to this puzzle are QCD's emer- 
gent phenomena: confinement and dynamical chiral symmetry break- 
ing. Confinement is the feature that no matter how hard one strikes a hadron, 
it never breaks apart into quarks and/or gluons that reach a detector. DCSB 
is signalled by the very unnatural pattern of bound state masses, something 
that is partly illustrated via Table 1 and the associated discussion. Neither of 
these phenomena is apparent in QCD's action and yet they are the dominant 
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determining characteristics of real- world QCD. Attaining an understanding of 
these phenomena is one of the greatest intellectual challenges in physics. 

In order to come to grips with DCSB it is first necessary to know the 
meaning of chiral symmetry. It is a fact that, at the Lagrangian level, local 
gauge theories with massless fermions possess chiral symmetry. Consider then 
hclicity, which may be viewed as the projection of an object's spin, j, onto 
its direction of motion, p ; viz., A oc j ■ p. For massless fermions, helicity is a 
Lorentz invariant spin observable. Plainly, it is cither parallel or anti-parallel 
to the direction of motion. 

In the Dirac basis, 75 is the chirality operator and one may represent a 
positive helicity (right-handed) fcrmion via 

q+(x) = \ (Id + 75) q{x) =: P+ q{x) (1.12) 
and a left-handed fermion through 

q-{x) = \ (Id - 75) q{x) =:P-q(x). (1.13) 

A global chiral transformation is enacted by 5 

q(x) -> q(x)' = c l ^ e q(x) , q(x) -> q(x)' = q{x) e 1 ^ 6 ', (1.14) 

and with the choice 9 = tt/2 it is evident that this transformation maps 
q+ — > g+ and g_ — > —q~. Hence, a theory that is invariant under chiral trans- 
formations can only contain interactions that are insensitive to a particle's 
hclicity; namely, helicity is conserved. 

Consider now a composite local pseudoscalar: q{x)i^^q{x). According to 
Eq. (1.14), a chiral rotation through an angle 9 = ir/4 effects the transforma- 
tion 

q~{x)i^ 5 q(x) -> -q(x)I D q{x) ; (1.15) 

i.e., it turns a pseudoscalar into a scalar. Thus the spectrum of a theory invari- 
ant under chiral transformations should exhibit degenerate parity doublets. 
Is such a prediction borne out in the hadron spectrum? Let's check (Rcf. [23], 
masses in MeV): 

tt(J p =0-,m= 140) cf. a(J p = 0+,m = 600) 

p(l~,TT0) cf. ai (l+,1260) (1.16) 

N(i + ,938) cf. A^(i",1535) 

Quite clearly, it is not: the difference in masses between parity partners is 
very large, which forces a conclusion that chiral symmetry is badly broken. 
Since the current-quark mass term is the only piece of the QCD Lagrangian 

5 For this illustrative purpose it is not necessary to consider complications that 
arise in connection with U(l) chiral anomalies, which appear via quantisation. 
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that breaks chiral symmetry, this appears to suggest that the quarks are quite 
massive. The conundrum reappears again: how can the pion be so light if the 
quarks are so heavy? 

The extraordinary phenomena of confinement and DCSB can be identified 
with properties of dressed-quark and -gluon propagators. These two-point 
functions describe the in-medium propagation characteristics of QCD's ele- 
mentary excitations. Here the medium is QCD's ground state; viz., the in- 
teracting vacuum. In propagating through this medium the quark and gluon 
propagators acquire momentum-dependent modifications that fundamentally 
alter the spectral properties of these elementary excitations. 

A mass term in the QCD Lagrangian explicitly breaks chiral symmetry. 
The effect can be discussed in terms of the quark propagator. For illustration 
it is sufficient to consider that of a nonintcracting fermion of mass m: 

S (P) = 2~T — — ■ ( L17 ) 

On this propagator, the chiral rotation of Eq. (1.14) is effected through 

S( P ) -> e** B S(p)e** B = Z 11 ' P 9 + e 2 ^ 9 9 m (1.18) 
p z + m 2 - p z + m l 

It is therefore clear that the symmetry violation is proportional to the current- 
quark mass and hence that the theory is chirally symmetric for m = 0. Another 
way of looking at this is to consider the fermion condensate: 

(99) = -tr / 7&S(p)oc - / p-^- 2 . (1.19) 
J (27r) 4 J (2np p A + m z 

This is a quantity that can rigorously be defined in quantum field theory [24] 
and whose strength measures the violation of chiral symmetry. It is a standard 
order parameter for chiral symmetry breaking, playing a role analogous to that 
of the magnetisation in a ferromagnet. 

This connects immediately to dynamical symmetry breaking. Consider a 
point-particle in the rotationally invariant potential 

V{<t,tt) = (cr 2 +7T 2 - l) 2 , (1.20) 

which is illustrated in Fig. 2. The figure depicts a state wherein the particle is 
stationary at an cxtremum of the action. That state is rotationally invariant 
but unstable. On the other hand, in the ground state of the system the particle 
is stationary at any point (cr, w) in the trough of the potential, for which a 2 + 
7r 2 = 1. There is an uncountable infinity of such vacua, \6), which are related, 
one to another, by rotations in the (a, 7r)-plane. The vacua are degenerate 
but not rotationally invariant and hence, in general, (0\cr\6) 7^ 7^ (8\ir\9). In 
this case the rotational invariance of the Hamiltonian is not exhibited in any 
single ground state: the symmetry is dynamically broken with interactions 
being responsible for (6\a\6) ^0^ (%|0). 
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Fig. 2. A rotationally invariant but unstable extremum of the Hamiltonian obtained 
with the potential in Eq. (1.20). (Adapted from Ref. [22].) 

The connection between dynamics and symmetries is now in plain view. 
The elementary excitations of QCD's action are absent from the strong inter- 
action spectrum: neither a quark nor a gluon ever reaches a detector alone. 
This is the physics of confinement. Chirality is the projection of a particle's 
spin onto its direction of motion. It is a Lorentz invariant for massless fermions. 
To classical QCD interactions, left-handed and right-handed quarks are indis- 
tinguishable. This symmetry has implications for the spectrum that do not 
appear to be realised. That is DCSB. Our challenge is to understand the emer- 
gence of confinement and DCSB from the QCD Lagrangian, and therefrom 
describe their impact on the strong interaction spectrum and hadron dynam- 
ics. These two phenomena need not be separate. They are both likely to be 
manifestations of the same mechanism. That mechanism must be elucidated. 
It is certainly nonperturbative. 

2 Nonperturbative Tool in the Continuum 

The Dyson-Schwinger equations (DSEs) provide a nonperturbative tool for 
the study of continuum strong QCD. At the simplest level they provide a 
generating tool for perturbation theory. Since QCD is asymptotically free, 
this means that any model-dependence in the application of these methods 
can be restricted to the infrared or, equivalently, to the long-range domain. 
In this mode, the DSEs provide a means by which to use nonperturbative 
strong interaction phenomena to map out, e.g., the behaviour at long range of 
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the interaction between light-quarks. A nonperturbative solution of the DSEs 
enables the study of: hadrons as composites of dressed-quarks and -gluons; the 
phenomena of confinement and DCSB; and therefrom an articulation of any 
connection between them. The solutions of the DSEs are Schwingcr functions 
and because all cross-sections can be constructed from such n-point functions 
the DSEs can be used to make predictions for real-world experiments. One 
of the merits in this is that any assumptions employed, or guesses made, 
can be tested, verified and improved, or rejected in favour of more promising 
alternatives. The modern application of these methods is described in Refs. [2, 
3, 4, 5, 6] and a pedagogical overview provided in Ref. 

Let's return to the dressed-quark propagator, which is given by the solution 
of QCD's gap equation that is expressed in Euclidean space as: 

S^ 1 =ij.p +m + S(p), (2.1) 

At zero temperature and chemical potential the most general Poincare co- 
variant solution of this gap equation involves two scalar functions. There are 
three common expressions: 

s{p) = viTmTW) = h-p+Mw = ~ n ' pav{p2) + asip2) ' 

(2.3) 

which arc equivalent. In the second form, Z{p 2 ) is called the wave-function 
renormalisation and M[p ) is the dressed-quark mass function. (A derivation 
of Eq. (2.1) and a discussion of its elementary and perturbative properties is 
presented in Sect. 2.2 of Ref. [5].) 

A weak coupling expansion of the DSEs produces every diagram in pertur- 
bation theory. The general result in perturbation theory can be summarised 



Sp 0r t (p ) 





\P 2 ' 






m 2 





(2.4) 



where the ellipsis denotes terms of higher order in a that involve (\n[p 2 / m 2 }) 2 
and (lnln[p 2 /m 2 ]), etc. However, at arbitrarily large finite order in perturba- 
tion theory it is always true that 

lim B port (p 2 ) =0. (2.5) 

m—*0 

This means that if one starts with a chirally symmetric theory, then in pertur- 
bation theory one also ends up with a chirally symmetric theory: the fermion 
DSE cannot generate a gap if there is no bare-mass seed in the first place. 
Thus DCSB is impossible in perturbation theory. Our question is whether this 
conclusion can ever be avoided; namely, are there circumstances under which 
it is possible to obtain a nonzero dressed-quark mass function in the chiral 
limit; viz., for m — > 0? 
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2.1 Dynamical Mass Generation 

To begin the search for an answer, consider Eqs. (2.1), (2.2) with the following 
model forms for the dressed-gluon propagator and quark-gluon vertex: 6 

g 2 D^{p-l) = 5 llv \e{A 2 -f), (2.6) 
m G 

r?(k,p)= 7 ^, (2.7) 

wherein m G is some gluon mass-scale and A serves as a cutoff. The model 
thus obtained is not renormalisable so that the regularisation scale A, upon 
which all calculated quantities depend, plays a dynamical role. In this case 
the gap equation is 

*7 ■ pA(p 2 ) + B(p 2 ) = i"f ■ p + mo 

4.J_ f dH 2 2 -* 7 -M(f 2 ) + g(l 2 ) ^ 

wherein mo is the mass that explicitly breaks chiral symmetry. 

If one multiplies Eq. (2.8) by (—i"f-p) and subsequently evaluates a matrix 
trace over spinor (Dirac) indices, then one finds 

2,81 f dH . A2 . 2 , „ A{P) 



P W = P + 3 m^ G J W ~ i)P ' e ¥AH^BW) • (2 ' 9) 

It is straightforward to show that / d 4 £p ■ I = 0; i.e., the angular integral in 
Eq. (2.9) vanishes, from which it follows that 

A(p 2 )^l. (2.10) 

This owes to the fact that models of the Nambu-Jona-Lasinio type are de- 
fined via a four-fermion contact-interaction in configuration space, which en- 
tails momentum-independence of the interaction and therefore also of the gap 
equation's solution in momentum space. 

If, on the other hand, one multiplies Eq. (2.8) by Jd, uses Eq. (2.10) and 
subsequently evaluates a trace over Dirac indices, then 



2^ 16 1 f d 4 

B{p 2 ) = m + — — — 
6 mi J (27T 



Since the integrand here is p 2 -independent then a solution at one value of p 2 
must be the solution at all values; viz., any nonzero solution must be of the 
form 



The form for the gluon two-point function implements a four-dimensional-cutoff 
version of the Nambu-Jona-Lasinio model, which has long been used to model 
QCD at low energies, e.g., Refs. [25, 26, 27]. This model stands within a class of 
four-fermion interaction models [28, 29]. 
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-BO 2 ) = constant = M . (2.12) 
Using this result, Eq. (2.11) becomes 

M = m +M-L \c{M 2 ,A 2 ), (2.13) 

07T 777 1 q 

C(M 2 , A 2 ) = A 2 - M 2 In [l + A 2 /M 2 ] . (2.14) 

Recall now that A defines the mass-scale in a nonrenormalisable model. Hence 
one can set A = 1 and hereafter merely interpret all other mass-scales as being 
expressed in units of A, whereupon the gap equation becomes 

M = m + M -L _L C(M 2 , 1) . (2.15) 
Let us consider Eq. (2.15) in the chiral limit: mo = 0, 

tf = 44 C(MV) - (2 ' 16) 

One solution is obviously M = 0. This is the result that connects smoothly 
with perturbation theory: one starts with no mass, and no mass is generated. 
In this instance the theory is said to realise chiral symmetry in the Wigner- 
Weyl mode. 

Suppose, however, that M ^ in Eq. (2.16). That is possible if, and only 
if, the following equation has a solution: 

1 = ^\C{M 2 ,1). (2.17) 

It is plain from Eq. (2.14) that C(M 2 , 1) is a monotonically decreasing function 
of M whose maximum value occurs at M = 0: C(0, 1) = 1. Consequently, 
3M 7^ solution if, and only if, 

i4 >x - (2 - l8) 

It is thus apparent that there is always a domain of values for the gluon mass- 
scale, rriQ, for which a nontrivial solution of the gap equation can be found. 
If one supposes A ~ 1 GeV, which is a scale that is typical of hadron physics, 
then 3M ^ solution for 



/12 

m 2 G < —r ~ (0.2 GeV) 2 . (2.19) 

07T 



This result, derived in a straightforward manner, is astonishing! It reveals 
the power of a nonperturbative solution to nonlinear equations. Although 
we started with a model of massless fermions, the interaction alone has pro- 
vided the fermions with mass. This is dynamical chiral symmetry breaking; 
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Fig. 3. With a bare mass m = 0.01, the mass gap obtained using the interaction 
specified by Eqs. (2.6), (2.7) - solid curve with circles. (All dimensioned quantities 
measured in units of A, the regularisation scale. Adapted from Ref. [5].) 



namely, the generation of mass from nothing. When this happens chiral sym- 
metry is said to be realised in the Nambu-Goldstone mode. It is clear from 
Eqs. (2.18), (2.19) that DCSB is guaranteed to be possible so long as the in- 
teraction exceeds a particular minimal strength. This behaviour is typical of 
gap equations. 

Figure 3 depicts the mc-dependence of the nontrivial self-consistent solu- 
tion of Eq. (2.15) obtained with a current-quark mass m = 0.01 (measured in 
units of A). The vertical line marks 



>G 



(2.20) 



namely, the value of mo above which a DCSB solution of Eq. (2.16) is impos- 
sible. Evidently, for ma > m G \ M(ma) ~ m and the self-consistent solution 
is well approximated by the perturbative result. However, a transition takes 
place for wiq ~ wg, and for mc < ttiq the dynamical mass is much greater 
than the bare mass, with M increasing rapidly as wiq is reduced and the 
effective strength of the interaction is thereby increased. 



2.2 Dynamical Mass and Confinement 



One aspect of quark confinement is the absence from the strong interaction 
spectrum of free-particle-like quarks. What does the model of Sect. 2.1 have 
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to contribute in this connection? Well, whether one works within the domain 
of the model on which DCSB takes place, or not, the quark propagator always 
has the form 

CNJL^ = 1 = -^-P + M 

yP> i 7 • p [A{p 2 ) = 1] + [B{p 2 ) =M] p 2 + M 2 ' 1 ' ' 

where M is constant. This expression has a pole at p 2 + M 2 = and thus is 
always effectively the propagator for a nonintcracting fcrmion with mass M. 
Hence, while it is generally true that models of the Nambu-Jona-Lasinio type 
support DCSB, they do not exhibit confinement. 

Consider an alternative [30] defined via Eq. (2.7) and 



g 2 D flu (k) = (2Tr) i GS i (k) 



On 



(2.22) 



Here G defines the model's mass-scale and the interaction is a (^-function in 
momentum space, which may be compared with models of the Nambu-Jona- 
Lasinio type wherein the interaction is instead a ^-function in configuration 
space. In this instance the gap equation is 

z 7 • p A(p 2 ) + B(p 2 ) = ll ■ p + m + (5 7m ^J^l^W ^ ' ^ 
which yields the following coupled nonlinear algebraic equations: 

a ^ =i+ vwW) ' (2 - 24) 

^ = ^ Vavw) ' (2 ' 25) 

Equation (2.23) yields an ultraviolet finite model and hence there is no rcgu- 
larisation mass-scale. In this instance one can therefore refer all dimensioned 
quantities to the model's mass-scale and set (5=1. 
Consider the chiral limit of Eq. (2.25): 

Obviously, like Eq. (2.16), this equation admits a trivial solution B(p 2 ) = 
that is smoothly connected to the perturbative result, but is there another? 
The existence of a B ^ solution; i.e., a solution that dynamically breaks 
chiral symmetry, requires (in units of G) 

p 2 A 2 (p 2 )+B 2 (p 2 )=A. (2.27) 



Suppose this identity to be satisfied, then its substitution into Eq. (2.24) gives 
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A{p 2 )-l= l -A{p 2 ) => A(p 2 )^2, (2.28) 



which in turn entails 



B{p 2 ) = 2^\-p 2 . (2.29) 

A complete chiral-limit solution is composed subject to the physical re- 
quirement that the quark self-energy is real on the spacelike momentum do- 
main, and hence 



2 ; p 2 < 1 



p 2 > 1 



(2.30) 



B(p 2 ) = {lf^'^4l (2-31) 

In this case both scalar functions characterising the dressed-quark propagator 
differ significantly from their free-particle forms and are momentum depen- 
dent. (We will see that this is also true in QCD.) It is noteworthy that the 
magnitude of the model's mass-scale plays no role in the appearance of this 
DCSB solution of the gap equation; viz., reinstating G, one has in the chiral 
limit B(0) = 2G ^ VG ^ 0. Thus, in models of the Munczek-Nemirovsky 
type, the interaction is always strong enough to support the generation of 
mass from nothing. 

The DCSB solution of Eqs. (2.30), (2.31) is defined and continuous for all 
p 2 , including timclike momenta, p 2 < 0. It gives a dressed-quark propagator 
whose denominator 

p 2 A 2 (p 2 ) + B 2 (p 2 ) >0, Vp 2 . (2.32) 

This is a novel and remarkable result, which means that the propagator 
does not exhibit any free-particle-like poles! It has long been argued that 
this feature can be interpreted as a realisation of quark confinement; e.g., 
Rcfs. [31, 32, 33, 34, 35, 36]. 

The Munczek-Nemirovsky interaction has taken a massless quark and 
turned it into something which at timelike momenta bears little resemblance 
to the perturbative quark. It does that for all nonzero values of the model's 
mass-scale. In this model one exemplifies an intriguing possibility that all 
models with quark-confinement necessarily exhibit DCSB. It is obvious from 
Sect. 2.1 that the converse is certainly not true. 

In the chirally asymmetric case the gap equation yields 

The second of these is a quartic equation for B(p 2 ). It can be solved alge- 
braically. There are four solutions, obtained in closed form, only one of which 
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Fig. 4. With a bare mass mo = 0.015, the dressed-quark mass function M(s = p 2 ) 
obtained using the interaction specified by Eqs. (2.7), (2.22) - solid curve. A free- 
particle-like mass-pole would occur at that value of a for which the solid and green 
curves intersect. As the figure suggests, this never happens in models of the Munczek- 
Nemirovsky type. (All dimensioned quantities measured in units of G, the model's 
mass-scale. Adapted from Ref. [5].) 

possesses the physically sensible ultraviolet spacelike behaviour: B(p 2 ) — > mo 
as p 2 — > oo. The physical solution is depicted in Fig. 4. At large spacelike mo- 
menta, M(s = p 2 ) — > mj and a perturbative analysis is reliable. That is 
never the case for s < 1 (in units of G), on which domain M(s) mo and 
the difference (M(s) — \s\) is always nonzero, a feature that is consistent with 
confinement. 

Two variations on the theme of dynamical mass generation via the gap 
equation have been illustrated. A general class of models for asymptotically 
free theories may be discussed in terms of an effective interaction 

g 2 D(Q 2 )=4n^±, (2.35) 

where a(Q 2 ) is such that g 2 D(Q 2 ) evolves according to QCD's renormalisation 
group in the ultraviolet. This situation is depicted in Fig. 5. The types of 
effective interaction fall into two classes. In those for which a(0) < 1, the only 
solution of the gap equation in the chiral limit is M(s) = 0. Whereas, when 
o(0) > 1, M(s) is possible in the chiral limit and, indeed, corresponds to 
the energetically favoured ground state. This behaviour was anticipated long 
ago in the context of attempts to understand QED at strong coupling [37, 38] . 
That theory unmodified, however, is incomplete [39]. 



Aspects of Hadron Physics 



17 



a(Q) 



0.2 - 



0.4 - 



0.8 



0.6 - 



1.0 - 




0.0 







10 



100 



Q (GeV) 



Fig. 5. Two classes of effective interaction in the gap equation. A theory will 
exhibit dynamical chiral symmetry breaking if, and only if, a(Q = 0) > 1. Here the 
ultraviolet behaviour of the interactions is constrained to be that of QCD. 

It is the right point to recall the conundrums described in Sect. 1. We have 
seen how confinement and DCSB can emerge in the nonperturbative solution 
of a theory's DSEs. Therefore theories whose classical Lagrangian possesses 
no mass-scale can, in fact, via DCSB, behave like theories with large quark 
masses and can also exhibit confinement. This understanding is a step toward 
a resolution of the riddles posed by strong interaction physics. 

2.3 Dynamical chiral symmetry breaking and a critical mass 

A novel aspect of the interplay between explicit and dynamical chiral symme- 
try breaking has recently been elucidated [40]. To relay this, suppose DCSB 
takes place so that B(p ) ^ 0. Then with the choice 8 = tt/2, Eq. (1.18) cor- 
responds to mapping B(p 2 ) — > —B{p 2 ). It follows that if B{p 2 ) is a solution 
of the gap equation in the chiral limit, then so is [— B(p 2 }. While these two 
solutions are distinct, the chiral symmetry entails that each yields the same 
pressure 7 and hence they correspond to equivalent vacua. This is an analogue 



7 The pressure is defined as the negative of the effective-action. Hence, the effective- 
action difference is zero between two vacuum configurations of equal pressure. A 
system's ground state is that configuration for which the pressure is a global max- 
imum or, equivalently, the effective-action is a global minimum. An elucidation 
may be found in Ref. [41]. 
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of the chiral-limit equivalence between the (a = 1, tt = 0) and [a = — 1,tt = 0) 
vacua in Fig. 2, as elucidated in Refs. [42, 43]. 8 

In common with Sects. 2.1 and 2.2, studies of DCSB have hitherto focused 
on a positive definite solution of the gap equation because the introduction of a 
positive current-quark bare-mass favours this solution; viz., if another solution 
exists, then it has a lower pressure. Returning again to Fig. 2, such a bare- 
mass tilts this so-called wine-bottle potential, producing a global minimum at 
(a = l,7r = 0). However, whether the massive gap equation admits solutions 
other than that which is positive definite, the effect of the current-quark mass 
on such solutions, and their interpretation, are questions little considered. 

Exemplar 

In order to begin addressing these questions, one may return to the model 
introduced in Sect. 2.1. The gap equation in this model, Eq. (2.15), can be 
rewritten as 



and, as noted above, anM^O solution exists when mo = if and only if vtlq 
satisfies Eq. (2.20). Hence, to proceed let us choose 



For this case G(M) is plotted in Fig. 6. One reads from the figure that in 
the chiral limit there are three solutions to Eq. (2.36): 



As observed after Eq. (2.16), Mw describes a realisation of chiral symmetry 
in the Wigncr mode. It corresponds to the vacuum configuration in which 
the possibility of DCSB is not realised; i.e., the position of the ball in Fig. 2. 
In the chiral limit this is the only solution accessible in perturbation theory. 
The solutions M± are essentially nonperturbative. Plainly, they represent the 
realisation of chiral symmetry in the Nambu-Goldstone mode; namely DCSB. 

M = M + > is the solution usually tracked in connection with QCD 
phenomenology; e.g., it is the solution plotted in Fig. 3. In models of this type 
M = 7U+ is identified as a constituent-quark mass. As mo is increased, M+ 
also increases. 



NB. More generally, given a solution of the mo > gap equation characterised 
by {^m (p 2 ),B mo (p 2 )}, then {A- mo (p 2 ), -B- mo (p 2 )} is a solution of the gap 
equation obtained with [—mo]. 




(2.36) 




(2.38) 
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Fig. 6. The zeros of G(M) give the solution of the gap equation defined by 
Eqs. (2.15), (2.37). Solid curve: obtained with mo = 0, in which case G(M) is odd 
under M — > —M; long-dashed curve: mo = 0.01; short-dashed curve: rno = m" = 
0.033; dotted curve: mo — 0.05. (All dimensioned quantities in units of A. Adapted 
from Ref. [40].) 

As evident in Fig. 6, the other two solutions of the gap equation do not im- 
mediately disappear when mo increases from zero. Nor do they always persist. 
Instead, these solutions exist on a domain 



Mw,M_ also evolve smoothly with m . Moreover, at the critical current- 
quark mass, m", these two solutions coalesce. 

To understand the origin of a critical mass one observes from Eq. (2.36) 
and Fig. 6 that introducing a current-quark mass merely produces a constant 
pointwise negative shift in the curve G{M). Hence, the critical current-quark 
mass is that value of this mass for which G(M) = at its local maximum. 
The local maximum occurs at 



and therefore in the present illustration the critical current-quark mass 



V(ma) = {mo | < mo < rriQ T }. 



(2.39) 



Mi 



{M | M < 0, G'{M)} = 



(2.40) 



m. 



cr 



— 



{m I G(M lm ) 



0} = 0.033 . 



(2.41) 



An interpretation of m" follows shortly. 
The two combinations 




(2.42) 
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Fig. 7. Evolution with current-quark mass of M (solid curve) and M (dashed 
curve). (All dimensioned quantities in units of A. Adapted from Ref. [40].) 



are of interest. In the chiral limit, M = M° and M = Mw- The evolution 
of each with current-quark mass is depicted in Fig. 7. It is apparent that M 
is continuous on 2?(mo) and evolves from the DCSB solution with increasing 
m Q : M(m ) is a monotonically decreasing function. These features can be 
understood as illustrating that the essentially dynamical component of chiral 
symmetry breaking decreases with increasing current-quark mass. 

The alternative combination. M, is also continuous on X>(mo). With the 
value of the coupling given in Eq. (2.37), M(mo) evolves from the Wigner 
solution according to 



M = mo 



1 



2 3(M ) 5 



1 



3 (M 



cm 



1 



(2.43) 



The development of M(mo) might be viewed as a gauge of the destabilising 
effect that DCSB has on this model's Wigner mode. 

What now can one make of the critical current-quark mass, Eq. (2.41)? 
Mw(mo) can be calculated perturbatively in the neighbourhood of mo = 0: 



M w {m Q ) = —3 mo + 



(2.44) 



However, with increasing current-quark mass Mw decreases steadily toward 
M-, which is an essentially nonperturbative quantity, until at m", Mw = 
M-. At this point a solution which began perturbatively has melded with a 
solution that is inaccessible in perturbation theory and actually characteristic 
of DCSB. A related view sees M_ as a DCSB solution whose modification by a 
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current-quark mass can no longer be evaluated pcrturbatively when that mass 
exceeds m". These observations suggest that to" specifies an upper bound 
on the domain within which, for physically relevant quantities, a perturbative 
expansion in the current-quark mass can be valid; i.e., it is a (possibly weak) 
upper bound on the radius of convergence for such an expansion. This view 
and the numerical result in Eq. (2.41) coincide with Ref. [44]. It is natural to 
ask whether analogous behaviour exists in QCD, and this is a question that 
will be addressed. First, however, additional concepts must be introduced. 



3 Gap equation in QCD 

The gap equation in QCD may properly be written 

S(p)- 1 = Z 2 {i7-p + m haIe )+Z 1 J g 2 D^(p-q)^S(q)r^q;p). (3.1) 

This is the renormalised DSE for the dressed-quark propagator wherein 
Dfj, v {k) is the renormalised dressed-gluon propagator, r^(q;p) is the renor- 
malised drcsscd-quark-gluon vertex, m bare is the vl-dcpendent current-quark 
bare mass that appears in the QCD Lagrangian and / := J d A q/(2ir) A rep- 
resents a Poincar e- invariant regularisation of the integral, with A the regular- 
isation mass-scale. In addition, Zi(C 2 , A 2 ) and Z 2 (( 2 , A 2 ) are the quark-gluon- 
vertex and quark wave function renormalisation constants, which depend on 
the renormalisation point, £, and the regularisation mass-scale. The solution 
of Eq. (3.1) is obtained subject to the renormalisation condition 

S(p) _1 | p2=c2 = iyp + m(Q, (3-2) 

where m(£) is the renormalised mass: 

Z 2 (C 2 , A 2 ) m bare (A) = Z A {C 2 , A 2 ) m(C) , (3.3) 

with Z/± the Lagrangian mass renormalisation constant. In QCD the chiral 
limit is unambiguously defined by 

Z 2 (C 2 , A 2 ) m haic (A) = , Vvl > C , (3.4) 

which is equivalent to stating that the renormalisation-point-invariant current- 
quark mass vanishes; i.e., to = 0. 

The discussion in Sect. 2 leads to the key question: what is the behaviour 
of the kernel of QCD's gap equation, Eq. (3.1). This kernel is constituted by 
the contraction of the dressed-gluon propagator and the dressed-quark-gluon 
vertex: 

D^{p-q)r y {q). (3.5) 
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Fig. 8. Solid curve - F(p 2 ) in Eq. (3.6) obtained through the solution of one 
particular truncation of the coupled ghost-gluon-quark DSEs with Nf = 3 flavours of 
massless quark. Dashed curve - DSE solution with the same truncation but omitting 
coupling to the quark DSE. Open circles - quenched lattice-QCD results. NB. There 
is sufficient lattice data at intermediate spacelike p 2 for the open circles to appear 
as a grey band. Within the lattice errors, the quenched DSE and lattice results are 
indistinguishable. (Adapted from Ref. [46].) 

In Landau gauge the two-point gluon Schwinger function can be expressed 



where F(p 2 ) involves the gluon vacuum polarisation discussed, e.g., in Sect. 2.1 
of Ref. [5]. The modern DSE perspective on F(p 2 ) is reviewed in Refs. [3, 6] and 
the predictions described therein were verified in contemporary simulations of 
lattice-regularised QCD [45]. The agreement is illustrated in Fig. 8. 

The DSE result depicted in Fig. 8 describes a gluon two-point function that 
is suppressed at small p 2 ; i.e., in the infrared. This deviation from expectations 
based on perturbation theory becomes apparent at p 2 ~ 1 GeV 2 . A mass-scale 
of this magnitude has long been anticipated as characteristic of nonperturba- 
tive gauge-sector dynamics. Rs origin is fundamentally the same as that of 
Aqqd, which appears in perturbation theory. This phenomenon whereby the 
value of a dimensionless quantity becomes essentially linked to a dynamically 
generated mass-scale is sometimes called dimensional transmutation. 

The dressed-quark two-point function has the form presented in Eq. (2.3). 
One knows that for a free particle Z(p 2 ) = 1 and M{p 2 ) = m currcnt . Asymp- 
totic freedom in QCD entails that at some large spacelike p 2 = £ 2 the quark 
propagator is exactly the free propagator except that the bare mass is replaced 




(3.6) 




0.0 1.0 2.0 3.0 4.0 0.0 1.0 2.0 3.0 4.0 

p (GeV) p (GeV) 

Fig. 9. Left Panel - Dashed-curve: gap equation's solution in the chiral limit; solid 
curves: solutions for M(p 2 ) obtained using the current-quark masses in Eq. (3.7). 
Data, upper three sets: lattice results for M(p 2 ) in GeV at am values in Eq. (3.7); 
lower points (boxes) : linear extrapolation of these results to am — [49] . Right Panel 
- Dashed curve, Z(p 2 ), and solid curve, M(p 2 ) calculated from the gap equation with 
m (C) = 55MeV. Data, quenched lattice-QCD results for M(p 2 ) and Z{p 2 ) obtained 
with am = 0.036 [49]. (Adapted from Ref. [51].) 

by the renormalised mass. At one-loop order, the vector part of the dressed 
self-energy is proportional to the running gauge parameter. In Landau gauge, 
that parameter is zero. Hence, in this gauge the vector part of the renor- 
malised dressed-quark self-energy vanishes at one-loop order in perturbation 
theory. The same is true for a charged fermion in QED. 

On the other hand, the behaviour of these functions on the nonperturbativc 
domain in QCD is a longstanding prediction of DSE studies [1] and could have 
been anticipated from Rcfs. [47, 48]. These DSE predictions, too, are confirmed 
in numerical simulations of lattice-QCD [49, 50], and the conditions have been 
explored under which pointwise agreement between DSE results and lattice 
simulations may be obtained [46, 51, 52]. This agreement is illustrated in 
Fig. 9, wherein the nonzero current-quark masses are 



a lattice 


0.018 0.036 0.072 


m(C) (GeV) 


0.030 0.055 0.110 



(3.7) 



The top row here lists the values used in the quenchcd-QCD lattice simu- 
lations, with a the lattice spacing so that am is dimcnsionless [49] and the 
second row provides the matched current-quark masses used in the DSE study, 
with the renormalisation scale C = Ci9 = 19 GeV [51]. 

Figure 9 confirms that DCSB is a reality in QCD. At ultraviolet momenta 
the magnitude of the mass function is determined by the current-quark mass. 
In the infrared, however, for light-quarks, M(p 2 ) is orders-of-magnitudc larger. 
The DSE analysis alone, and its correlation of lattice data, indicates that the 
mass function is nonzero and retains its magnitude in the chiral limit. In 
Fig. 9 and subsequently Fig. 12 the evolution in QCD from current-quark to 
constituent-quark is in plain view. 
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What about confinement? It is notable that confinement between static 
quarks is a quantum mechanical archetype. Numerical simulations of lattice- 
regularised QCD exhibit a linear potential at long range between such in- 
finitely heavy sources/sinks of colour charge. The physics of confinement in 
the light-quark domain, however, is not obviously amenable to such a de- 
scription. For example, a veracious realisation of dynamical chiral symmetry 
breaking is impossible in quantum mechanical potential models. Moreover, 
in the presence of light- (dynamical-) quarks, the breaking of the string that 
characterises the potential between static quarks appears to be an instan- 
taneous process, with delocalised light-quark-antiquark pair production [53]. 
Thus either no information is to be had about a string between light-quarks, 
or such a string cannot exist. 

Another perspective on confinement has already been mentioned; namely, 
that this phenomenon might be expressed in the analyticity properties of the 
dressed propagators. In fact it is sufficient for confinement that dressed prop- 
agators for coloured excitations not possess a Lehmann representation, since 
this is associated with a violation of reflection positivity. An excitation con- 
nected with a propagator that violates reflection positivity cannot appear in 
the Hilbcrt space of physical states; viz., it won't propagate to a detector. 
These notions may be traced from Refs. [31, 32, 33, 34, 35, 36] and are de- 
scribed in Refs. [1, 2, 3]. (See also, e.g., the discussion of the reconstruction 
theorem in Ref. [54].) Any Schwinger function that exhibits an inflexion point 
cannot be expressed through a Lehmann representation. An inspection of the 
DSE and lattice results for F(p 2 )/p 2 and <7s{p 2 ) suggests strongly that the 
dressed- gluon and -quark propagators display an inflexion point. Moreover, in 
QCD there are DSE studies which suggest that quark confinement and DCSB 
both owe to the same dynamical mechanism [55, 56], and therefore that one 
does not appear without the other. 

In connection with Fig. 5 in Sect. 2.2 a description was provided of a 
class of effective interactions in the gap equation that can generate DCSB. 
The dressed-gluon propagator obtained from F(p 2 ) in Fig. 8, combined with 
Eq. (2.7), does not yield an interaction that is a member of that class 
[57, 58, 59]. How then is it possible that Ref. [51] unified the gluon and quark 
two-point functions? The problem was anticipated in Ref. [59] and an answer 
suggested; namely, the dressed-quark-gluon vertex must exhibit an enhance- 
ment in the infrared. This is precisely the means employed in Ref. [51]. The 
exact nature of the enhancement and its origin in QCD is the subject of 
contemporary research [52, 60, 61]. 

4 Bound States 

4.1 Bethe-Salpeter Equation 

It will now be apparent that a semi-quantitatively reliable picture of QCD's 
key propagators and vertices is established. What about bound states? With- 
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out them, of course, a direct comparison with experiment is impossible. Bound 
states appear as pole contributions to colour-singlet Schwingcr functions and 
this observation may be viewed as the origin of the Bethe-Salpeter equation. 
The Bethe-Salpeter equation (BSE) has an history that predates QCD, which 
can be traced, e.g., from Rcf. [9], Chap. 10. 

The axial- vector vertex is of primary interest to hadron physics. It may be 
obtained as the solution of the inhomogeneous Bethe-Salpeter equation 

[/V(fc;P)] tu = Z 2 [ 757M ] tu + [ A {S(q + )r 5tl ( q ;P)S( q ^} sr K™(q,k;P) , (4.1) 

Jq 

where q± = q ± P/2 and the colour-, Dirac- and flavour-matrix structure of 
the elements in the equation is denoted by the indices r,s,t,u. In Eq. (4.1), 
K (g, k; P) is the fully-amputated quark-antiquark scattering kernel. It is one- 
particle-irreducible and hence, by definition, does not contain quark-antiquark 
to single gauge-boson annihilation diagrams, such as would describe the lep- 
tonic decay of the pion, nor diagrams that become disconnected by cutting one 
quark and one antiquark line. If one knows the form of K then one completely 
understands the nature of the interaction between quarks in QCD. 

4.2 Model-independent results 
Pseudoscalar meson mass formula 

Much has been made of chiral symmetry in the preceding discussion. In quan- 
tum field theory, chiral symmetry and the pattern by which it is broken are 
expressed via the chiral Ward-Takahashi identity: 

p^rilik: P) = s(k + )-H l5 ^- + i75 ^^ (fc _)-i _ i {M c j F H [k . P)}) (42) 

where the pseudoscalar vertex is given by 

+ I* [ X f(q;P)] sr KZ(q,k;P), (4.3) 
tu J q 

with xf(<?;P) = S(q+)ri I (q;P)S(q-), S = diag[5 u , S d , S s , . . .} and M< - 
diag[m u (C),m d (C),m s (C), . . .]. 

Equations (4.2), (4.3) are written for the case of a flavour-nonsinglct vertex 
in a theory with Nf quark flavours. The matrices T H are constructed from the 
generators of SU(Nf) with, e.g., T 7r+ = |(A X + iX 2 ) providing for the flavour 
content of a positively charged pion. Writing the equations in this manner 
is straightforward. However, a unified description of light- and heavy-quark 
systems is not. Truncations and approximations that are reliable in one sector 
need not be valid in the other. 



[ r 5 H (k;P)] tu = Z 4 



75" 
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The axial-vector Ward-Takahashi identity relates the solution of a BSE 
to that of the gap equation. If the identity is always to be satisfied and in 
a model-independent manner, as it must be in order to preserve an essential 
symmetry of the strong interaction and its breaking pattern, then the kernels 
of the gap and Bethe-Salpeter equations must be intimately related. Any 
truncation or approximation of these equations must preserve that relation. 
This is an extremely tight constraint. Perturbation theory is one systematic 
truncation that, order by order, guarantees Eq. (4.2). However, as has been 
emphasised, perturbation theory is inadequate in the face of QCD's emergent 
phenomena. Something else is needed. 

Fortunately, at least one systematic, nonperturbative and symmetry pre- 
serving truncation of the DSEs exists [61, 62, 63]. This makes it possible to 
prove Goldstone's theorem in QCD [64]. Namely, when chiral symmetry is 
dynamically broken: the axial- vector vertex, Eq. (4.1), is dominated by the 
pion pole for P 2 ~ and the homogeneous, isovector, pseudoscalar BSE has 
a masslcss (P 2 — 0) solution. The converse is also true, so that DCSB is 
a sufficient and necessary condition for the appearance of a masslcss pseu- 
doscalar bound state of dynamically-massive constituents which dominates 
the axial- vector vertex for infrared total momenta. 

Furthermore, from the axial-vector Ward-Takahashi identity and the exis- 
tence of a systematic, nonperturbative symmetry-preserving truncation, one 
can prove the following identity involving the mass-squared of a pseudoscalar 
meson [64]: 

f H m 2 H = p H (0M c H , (4.4) 

where Mjj = m qi (£) + m q2 (£) is the sum of the current-quark masses of the 
meson's constituents; 

f H P^ = Z 2 trj ^T H ) T 75 ^S(q + )r H (q;P)S(q-), (4.5) 

where (-) T indicates matrix transpose, and 

p H (0 = Z i trj* 1 - (T H f l5 S(U )r H (q;P)S(q_ ). (4.6) 

The renormalisation constants in Eqs. (4.5), (4.6) play a pivotal role. Indeed, 
the expressions would be meaningless without them. They serve to guarantee 
that the quantities described are gauge invariant and finite as the regulari- 
sation scale is removed to infinity, which is the final step in any calculation. 
Moreover, Z 2 in Eq. (4.5) and Z4 in Eq. (4.6) ensure that both fn and the 
product pniQMjj are renormalisation point independent, which is an abso- 
lute necessity for any observable quantity. 

Taking note that in a Poincare invariant theory a pseudoscalar meson 
Bethe-Salpeter amplitude assumes the form 
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ir J H {k; P) = T H l5 [iE H (k; P) + 7 ■ PF H (k; P) 
+ 1 -kk-PG H {k-P)+a^ k^ H H (k; P)} , 



(4.7) 



then, in the chiral limit, one can also prove that 



F R (k;0) + 2f H F H (k;Q) 
G R (k;0) + 2f H G H (k;0) 
H R (k;0) + 2f H H H (k;Q) 



f H E H (k;0) 



B(k 2 ), 
A(k 2 ), 
2A'(k 2 ) , 
0. 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



The functions F R , G R , H r arc associated with terms in the axial- vector vertex 
that are regular in the neighbourhood of P 2 + m 2 H = and do not vanish at 
Pft = 0. These four identities are quark-level Goldbcrgcr-Treiman relations for 
pscudoscalar mesons. They are exact in QCD and are a pointwise expression 
of Goldstone's theorem. These identities relate the pseudoscalar meson Bethe- 
Salpcter amplitude directly to the dresscd-quark propagator. Equation (4.8) 
explains why DCSB and the appearance of a Goldstone mode are so intimately 
connected, and Eqs. (4.9) - (4.11) entail that in general a pseudoscalar meson 
Bcthc-Salpeter amplitude has what might be called pseudovector components; 
namely: Fr, Gr, Hr- The latter are connected with the presence of quark 
orbital angular momentum in the pion, which is a necessary feature of a 
Poincare covariant treatment, and guarantee that the electromagnetic pion 
form factor behaves as 1/Q 2 at large spacelike momentum transfer [65, 66]. 

Equation (4.4) and its corollaries arc of fundamental importance in QCD. 
To exemplify let's focus first on the chiral limit behaviour of Eq. (4.6) whereat, 
using Eqs. (4.7), (4.8)-(4.11), one finds readily 



where is the chiral limit value from Eq. (4.5), which is nonzero when chiral 
symmetry is dynamically broken. Equation (4.12) is unique as the expression 
for the chiral limit vacuum quark condensate, and is the true definition of the 
order parameter first described in Eq. (1.19). It thus follows from Eqs. (4.4), 
(4.12) that in the neighbourhood of the chiral limit 



Hence what is commonly known as the Gell-Mann-Oakes-Renner relation is 
a corollary of Eq. (4.4). 

Let's now consider another extreme; viz., when one of the constituents 
is a heavy quark, a domain on which Eq. (4.4) is equally valid. In this case 
Eq. (4.5) yields the model-independent result [67] 





(4.12) 




M C H (qq}° ( + 0(M 2 ) . 



(4.13) 



1 



(4.14) 
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i.e., it reproduces a well-known consequence of heavy-quark symmetry [68]. A 
similar analysis of Eq. (4.6) gives a new result [69, 70] 



where ttiq is the renormalisation-group-invariant current-quark mass of the 
flavour-nonsinglet pseudoscalar meson's heaviest constituent. This is the re- 
sult one would have anticipated from constituent-quark models but here is 
indicated a direct proof in QCD. 

Pseudoscalar meson excited states 

Pseudoscalar mesons hold a special place in QCD and there are three states, 
composed of u,d quarks, in the hadron spectrum with masses below 2 GcV 
[23]: 7r(140); tt(1300); and tt(1800). Of these, the pion [tt(140)] is naturally 
well known and much studied. The other two are observed, e.g., as resonances 
in the coherent production of three pion final states via pion-nucleus collisions 
[71]. In the context of a model constituent-quark Hamiltonian, these mesons 
are often viewed as the first three members of a QQ n ^So trajectory, where n is 
the principal quantum number; i.e., the 7r(140) is viewed as the S'-wave ground 
state and the others are its first two radial excitations. By this reasoning the 
properties of the 7r(1300) and 7r(1800) are likely to be sensitive to details of the 
long-range part of the quark-quark interaction because the constituent-quark 
wave functions will possess material support at large intcrquark separation. 9 
Hence the development of an understanding of their properties may provide 
information about light-quark confinement, which complements that obtained 
via angular momentum excitations [73, 74]. As has been emphasised, the de- 
velopment of an understanding of confinement is one of the most important 
goals in physics. 

We have already seen that Eq. (4.4) is a powerful result. That is further 
emphasised by the fact that it is applicable here, too [75, 76]. The result 
holds at each pole common to the pseudoscalar and axial-vector vertices and 
therefore it also impacts upon the properties of non ground state pseudoscalar 

9 One might note that, in comparison with 7r(1300), the 7r(1800) is narrow, with a 
width of 207 ± 13 MeV [23], and has a decay pattern that may be consistent with 
its interpretation as a hybrid meson in constituent-quark models [72]. This picture 
has the constituent-quarks' spins aligned to produce Sqq = 1 with J — ob- 
tained by coupling Sqq to a spin-1 excitation of the confinement potential. From 
this perspective, too, the properties of 7r(1800) are sensitive to the interaction 
responsible for light-quark confinement. 




(4.15) 



Combining Eqs. (4.14), (4.15), one finds [69, 70] 



win oc rhn for 



(4.16) 
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mesons. Let's work with a label n > for the pseudoscalar mesons: 7r n , with 
n = denoting the ground state, n = 1 the state with the next lowest mass, 
and so on. Plainly, by assumption, Ti% nn . > m I0 , and hence m^ n ^ a > in 
the chiral limit. Moreover, the ultraviolet behaviour of the quark- antiquark 
scattering kernel in QCD guarantees that Eq. (4.6) is cutoff independent. Thus 

< p° n (C) := lim p w „(C) < oo, Vn. (4.17) 

m— >0 

Hence, it is a necessary consequence of chiral symmetry and its dynamical 
breaking in QCD; viz., Eq. (4.4), that 

# n = 0,Vn>l. (4.18) 

This result is consistent with Refs. [77] , as appreciated in Rcf. [78] . It means 
that in the presence of DCSB all pseudoscalar mesons except the ground state 
decouple from the weak interaction. 

This argument is legitimate in any theory that has a valid chiral limit. It is 
logically possible that such a theory does not exhibit DCSB; i.e., realises chiral 
symmetry in the Wigner-Weyl mode, as was illustrated in Sect. 2.1 and is the 
case in QCD above the critical temperature for chiral symmetry restoration 
[2]. Equation (4.4) is still valid in the Wigner phase. However, its implications 
are different; namely, in the Wigner phase, one has 

B w (0,C 2 ) ocm(C) ocm; (4.19) 

i.e., the mass function and constituent-quark mass vanish in the chiral limit. 
This result is accessible in perturbation theory. Equation (4.8) applies if there 
is a masslcss bound state in the chiral limit. Suppose such a bound state 
persists in the absence of DCSB. It then follows from Eqs. (4.8) and (4.19) 
that 

(4-20) 

In this case the leptonic decay constant of the ground state pseudoscalar 
meson also vanishes in the chiral limit, and hence all pseudoscalar mesons are 
blind to the weak interaction. 
It is always true that 

U Pn (0^°-(m)l (4.21) 

In the Wigner phase [24], (qq)® w <x to 3 . Hence, via Eq. (4.4), if a rigor- 
ously chirally symmetric theory possesses a masslcss pseudoscalar bound state 
then 10 

t° to . (4.22) 

In this case there is also a degenerate scalar meson partner whose mass behaves 
in the same manner. 



10 cf. The case of DCSB; i.e., the Nambu phase, wherein rn^ o "oc° y/rh, Eq. (4.13). 
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Pseudoscalar meson two photon decays 

In the presence of DCSB the ground state neutral pseudoscalar meson de- 
cays predominantly into two photons, a process connected with the Abelian 
anomaly. Given that /tt^o = in the chiral limit, it is natural to ask whether 
7r„^o — > 77 is similarly affected. Since rainbow-ladder is the leading order 
in a symmetry preserving truncation, it can be used to provide a model- 
independent analysis of this process. In this instance, the two-photon cou- 
pling for all u,d pseudoscalar mesons, including n = 0, is described by the 
renormalised triangle diagrams 

i-M 

j i 

x iQr^£-o,£)S{e)iQr u (i,e 0+ ), (4.23) 

i-M 

Tl^iM = tr J S(to + )r3(e Q+ ,£-o)S(£-o) 

x iQr^t-oJ) S(£)iQr v (£,£ 0+ ), (4.24) 

where £ a p = £ + ak\ + (3k2, the electric charge matrix Q = diag[e u ,ed] = 
ediag[2/3, -1/3], S = di&g[S u , S d ] and 

[r„(fc; P)] tu = Z 2 [ ltl ] tu + [ [£(g+)r„(g; P)$( q _)} sr K%( q , k; P) (4.25) 

Jq 

is the renormalised dressed-quark-photon vertex. 11 Is is subsequently assumed 
that m u (C.) = m(Q = m,j(C) so that S u = Sd- 

The dressed-quark propagators in Eqs. (4.23) - (4.25) are understood to 
be calculated using a rainbow-truncation gap equation, and if is a ladder-like 
quark-antiquark scattering kernel that ensures Eq. (4.2) is satisfied. At this 
point , nothing more need be said of these elements. The results described are 
independent of the details discussed in Sect. (4.3). 

Under the conditions just described one may derive [80] the Ward-Takahashi 
identity for the quark-antiquark scattering matrix, G, depicted in Fig. 10. Us- 
ing this identity it can be shown [81] 

P p Tl lvp (ki,k 2 ) + 2zm(C) Tl^M) = £ ^ vpa k lp k 2<7 , (4.26) 

where a = e 2 /(Air). This is an explicit demonstration that the triangle- 
diagram representation of the axial-vector-two-photon coupling calculated in 
the rainbow-ladder truncation is a necessary and sufficient pairing to preserve 
the Abelian anomaly. 

11 This is an inhomogeneous BSE for the dressed-quark-photon vertex. Its lowest 
mass pole-contribution is the ground-state p-meson [79]. That fact underlies the 
success of vector meson dominance phenomenology. 
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Fig. 10. This axial- vector Ward-Takahashi identity is an analogue of Eq. (4.2). 
It is valid if, and only if: the dressed-quark propagator, S, is obtained from 
Eq. (4.46); the axial- vector vertex, i"5 M , is obtained from Eq. (4.1) with a ladder- 
like kernel, K, that ensures Eq. (4.2); the pseudoscalar vertex is constructed anal- 
ogously; and the unamputated renormalised quark-antiquark scattering matrix: 
G = (SS) + {SS)K(SS) + (SS)K(SS)K(SS) +[...], is constructed from the el- 
ements just described. (Adapted from Ref. [76].) 



In general the coupling of an axial-vector current to two photons is de- 
scribed by a six-point Schwinger function, to which Eq. (4.23) is an approx- 
imation. The same is true of the pseudoscalar-two-photon coupling and its 
connection with Eq. (4.24). Equation (4.26) is valid for any and all values of 
P 2 = (fci + fc2) 2 - It is an exact statement of a divergence relation between these 
two six-point Schwinger functions, which is preserved by the rainbow-ladder 
truncation, and any systematic improvement thereof. It is now possible to 
report corollaries of Eq. (4.26) that have important implications for the prop- 
erties of pseudoscalar bound states. 

It follows from the discussion presented above that, unless there is a reason 
for the residue to vanish, every isovector pseudoscalar meson appears as a pole 
contribution to the axial- vector and pseudoscalar vertices [64] : 



P 2 +771,2 p2 



m 



iri(k;P) = f" (C) 2 rl(k;P)+ irt*(k;P), (4.28) 

where the residues are given in Eqs. (4.5), (4.6); viz., each vertex may be 
expressed as a simple pole plus terms regular in the neighbourhood of this 
pole, with ri (k; P) being the bound state's canonically normalised Bethe- 
Salpcter amplitude. 12 

If one inserts Eqs. (4.27) and (4.28) into Eq. (4.26), and uses Eq. (4.4), it 
follows that in the neighbourhood of each electric-charge-neutral pscudoscalar- 
meson bound-state pole 



P P T 5^( k i> + 2m(C) T^° s (k 1: k 2 ) + U n T^(k 1: k 2 ) = —ie t _ lly p a k lp k 2 a ■ 

(4.29) 

12 Canonical normalisation of the bound-state Bethe-Salpeter amplitude is dis- 
cussed, e.g., in Sect. 10-2-1 of Ref. [9] and Ref. [82]. The normalisation condition 
assumes a simple form when the two-body scattering amplitude is independent 
of the total momentum, as illustrated in Eq. (5.34). 
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In this equation, T 3reg (/ci, k 2 ) are nonrcsonant or continuum contributions to 
the relevant Schwingcr functions, whose form is concretely illustrated upon 
substitution of r^ cs (k;P) and r| rog (fc;P) into Eqs. (4.23) and (4.24), re- 

spectively. Moreover, T M! ? is the six-point Schwingcr function describing the 
bound state contribution, which in rainbow-ladder truncation is realised as 

rM-HX> 

3]£?(fci,fa)=tr / S(i 0+ )r < (i_ ii ;P)S(£-o) 

x iQr fi (t_ ,e)g(£)iQr„(£,e 0+ ). (4.30) 

This Schwingcr function describes the direct coupling of a pseudoscalar meson 
to two photons. The support properties of the bound state Bethc-Salpeter 
amplitude guarantee that the renormalised Schwingcr function is finite so 
that the regularising parameter can be removed; i.e., M — ► oo. 

Owing to the 0(4) (Euclidean Lorentz) transformation properties of each 
term on the l.h.s. in Eq. (4.26), one may write 



PpT^t B P (ki,k 2 ) = -ie pypa k lp k 2a A 3 ^(h,k 2 ) , (4.31) 
T tt S (h, ha) = -ie pvp<y k lp k 2<7 P^{kr,k 2 ) , (4.32) 

^ 7T 

Q/ o 

Tj^fci, fc 2 ) = -ie pl/pa k lp k 2a G^^iki, k 2 ) , (4.33) 

so that Eq. (4.26) can be compactly expressed as 

A 3 "*(fci,fca) + 2im(C)P 3res (A;i,fe) + U n G<(k u k 2 ) = i (4.34) 

It follows from Eqs. (4.18), (4.29) that in the chiral limit all pseudoscalar 
mesons, except the Goldstone mode, decouple from the divergence of the axial- 
vector-two-photon vertex. 

In the chiral limit the pole associated with the ground state pion appears 
at P 2 = and thus 



p P T 5^p(ki, fe)| F2 ^ = P P T 5 ™ g p (ki,k 2 ) = — ie pvpa k lp k 2<7 ; (4.35) 

namely, outside the neighbourhood of the ground state pole the regular (or 
continuum) part of the divergence of the axial-vector vertex saturates the 
anomaly in the divergence of the axial-vector-two-photon coupling. On the 
other hand, in the neighbourhood of P 2 = 

A 3 ™z(h,k 2 )\ p2 ^ + U G*°(k u k 2 ) = i; (4.36) 

i.e., on this domain the contribution to the axial- vector-two-photon coupling 
from the regular part of the divergence of the axial-vector vertex combines 
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with the direct 7Tq77 vertex to fulfill the anomaly. This fact was illustrated in 
Ref. [66] by direct calculation: Eqs. (4.8) - (4.11) are an essential part of that 
demonstration. 
If one defines 



T>(P 2 ,Q 2 ) = G^{kiM) 



k \ — Q ^ — k 



(4.37) 



in which case P 2 = 2(ki ■ k% + Q 2 ), then the physical width of the neutral 
ground state pion is determined by 



m w o,0); 



(4.38) 



viz., the second term on the l.h.s. of Eq. (4.36) evaluated at the on-shcll points. 
This result is not useful unless one has a means of estimating the contribu- 
tion from the first term; viz., A 3reg (fc 1 , k 2 ). However, that is readily done. A 
consideration [64] of the structure of the regular piece in Eq. (4.27) indicates 
that the impact of this continuum term on the 7Tq77 coupling is modulated by 
the magnitude of the pion's mass, which is small for realistic u and d current- 
quark masses and vanishes in the chiral limit. One therefore expects this term 
to contribute very little and anticipates from Eq. (4.36) that 



1 1 

2 /tt 



(4.39) 



is a good approximation. This is verified in explicit calculations; e.g., in 
Ref. [83] , which evaluates the triangle diagrams described herein, the first 
term on the l.h.s. modifies the result in Eq. (4.39) by less than 2%. 

There is no reason to expect an analogous result for pseudoscalar mesons 
other than the 7r(140); i.e., the states which are denoted by n > 1. Indeed, as 
all known such pseudoscalar mesons have experimentally determined masses 
that are greater than 1 GcV, the reasoning used above suggests that the pres- 
ence of the continuum terms, A 3Tes (ki, k^) and P 3reg (A;i, k%), must materially 
impact upon the value of <7 7r o 77 - This was shown to be true in Ref. [76]; e.g., 

It was stated that the rainbow-ladder truncation preserves the one-loop 
renormalisation group properties of QCD. It follows that Eq. (4.30) should 
reproduce the leading large-Q 2 behaviour of the J* (Q)^n(P)j* (Q) transition 
form factor inferred from perturbative QCD. That is true [84]. 

Reference [76] considered the general case; i.e., arbitrary n, and proved 



2 Q 2 »j4qcd 4tt 2 



Q 2 " 1 



Q 4 



(4.40) 



where: 7 is an anomalous dimension, which cannot be determined accurately 
in rainbow-ladder truncation; w, n is a mass-scale that gauges the momen- 
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turn space width of the pseudoscalar meson; and Fn (— Tn% ) is a structure- 
dependent constant, similar but unrelated to f nn . 13 It is now plain that V n > 1 



lun T^o (-m Q ) = — F 7 \ > (-m ) — a 

m— >0 o C,/ 



; (4.41) 

m=0 



namely, in the chiral limit the leading-order power-law in the transition form 
factor for excited state pseudoscalar mesons is 0(1/Q 4 ). This is a novel, non- 
perturbative and model- independent result. 

Reference [76] also provides quantitative illustrations of these and other 
results relating to the electromagnetic properties of ground- and excited-state 
pseudoscalar mesons. For example, the two photon width of the 7Ti-meson is 
estimated 

T X177 ~ 240 cV ~ 30 , (4-42) 

as is the electromagnetic radius of the charge states 

rv, ~ 1.4rv . (4.43) 

These results were obtained with the model rainbow-ladder interaction 
described below, in Sect. 4.3. However, the calculation of the latter was used 
to verify that the properties of excited states are indeed a sensitive probe of 
the interaction between light-quarks at long-range. This is shown in Fig. 11. 
As with so many properties of the ground state Goldstone mode, its radius 
is protected against rapid evolution by DCSB. However, the charge radius 
of the first excited state changes swiftly with increasing ui, with the ratio 
fwi/fwo varying from 0.9 - 1.2. This outcome can readily be interpreted. The 
length-scale r a := 1/lu is a measure of the range of strong attraction: magni- 
fying r a increases the active range of the confining piece of the interaction. 
This effectively strengthens the confinement force and that compresses the 
bound state, as one observes in Fig. 11: r ni decreases quickly with increasing 
r a = 1/uj. While it is natural to suppose r ni > r^ ; namely, that a radial exci- 
tation is larger than the associated ground state, the calculations of Ref. [76] 
illustrated that with the ground state pseudoscalar meson's properties con- 
strained by Goldstone's theorem and its pointwise consequences, Eqs. (4.8) 
- (4.11), it is possible in quantum field theory for a confining interaction to 
compress the excited state with the consequence that r 7ri < r Wo . 

A related result for the evolution of the mass was observed in Ref. [75]; 
namely, the mass of the first excited state dropped rapidly with increasing 
r a . (However, in this case DCSB guarantees m Tl > m To .) On the w-domain 
illustrated in Fig. 11, the mass of the ground state obtained with nonzero 
current-quark mass varied by only 3% while that of the first excited state 
changed by 14%. It is natural to expect that an increase in the strength of the 
confinement force should increase the magnitude of the binding energy and 
hence reduce the mass. That is precisely what occurs. 



13 With the interaction of Eq. (4.44), F^\-m% n ) ~ -(qq)°, and it is generally 
nonzero in the chiral limit. 
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Fig. 11. Evolution of ground and first excited state pseudoscalar mesons' electro- 
magnetic charge radii with the scale parameter ui in Eq. (4.44): r a — 1/w gauges the 
range of the confining interaction between light quarks. Dotted line: r„ = 0.66 fm, 
which indicates the experimental value of the ground state's radius. The radius is 
evaluated numerically from the electromagnetic form factor. That is the origin of 
the theory error depicted in the figure, which corresponds to a relative error < 1% 
for n = and < 3% for n = 1. (Adapted from Ref. [76].) 



4.3 Practical, predictive tool 

It was stated that there is a systematic, nonperturbative symmetry-preserving 
truncation of the DSEs. The leading-order is provided by the renormalisation- 
group-improved rainbow-ladder truncation, which has been used widely; e.g., 
Refs. [85, 86, 87], and references thereto. A practical renormalisation-group- 
improved rainbow-ladder truncation preserves the one-loop ultraviolet be- 
haviour of perturbative QCD. However, a model assumption is required for 
the behaviour of the kernel in the infrared; viz., on the domain Q 2 < 1 GcV 2 , 
which corresponds to length-scales > 0.2 fm. This is the confinement domain 
whereupon little is truly known about the interaction between light-quarks. 
That information is, after all, what we seek! 

The systematic application of a single model to an extensive range of JLab- 
related phenomena is pursued in Refs. [79, 83, 88, 89, 90, 91, 92, 93, 94]. The 
model interaction is expressed via 

^ = Q 2 e"^ 2 + — v ^ F(Q 2 ) , (4.44) 



1 + 2 M 2 qcd" 
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wherein F(Q 2 ) = [1 - exp(-Q 2 /[4m 2 ])]/Q 2 , m t = 0.5 GcV; r = c 2 - 1; 7m = 
12/25; and A QC d = = 0.234 GeV. (NB. Eq. (4.44) gives a(m|) = 0.126.) 
This simple form represents the interaction strength as a sum of two terms. 
The second ensures that perturbative behaviour is preserved at short-range. 
The first, on the other hand, provides for the possibility of enhancement in 
K and the gap equation's kernel at long-range [51]. The model assumption 
explicit here is restricted to the infrared; viz., to the interval Q 2 < 1 GeV 2 , 
which corresponds to length-scales > 0.2 fm. 

The true parameters in Eq. (4.44) are D and co, which together determine 
the integrated infrared strength of the rainbow-ladder kernel; i.e., the so-called 
interaction tension o A [58]. However, it is worth emphasising that they are 
not independent [88]: in fitting to a selection of observables, a change in one 
is compensated by altering the other; e.g., on the domain ui € [0.3,0.5] GeV, 
the fitted observables are approximately constant along the trajectory [95] 

u) D = (0.72 GcV) 3 =: m 3 . (4.45) 

This correlation: a reduction in D compensating an increase in lo, acts to 
keep a fixed value of the interaction tension. Equation (4.44) is thus a one- 
parameter model. NB. As we saw in connection with Fig. 8, this value of m g 
is typical of the mass-scale associated with nonperturbative gluon dynamics. 

With the interaction in Eq. (4.44), one obtains the rainbow-truncation gap 
equation: 

Sip)- 1 = Z 2 (z TP +m h n+J q A G((p-<l) 2 )D l * c u c (p-q) Y^S(q) y>, (4.46) 

wherein D l ™°(£) is the free gauge boson propagator; namely, Eq. (3.6) with 
F(p 2 ) = 1. The self-consistent solution obtained with current-quark masses 
(in GcV): 

(4.47) 

is depicted in Fig. 12. It is important to bear in mind and evident in the figure 
that, because the truncation preserves the one- loop renormalisation group 
properties of QCD, the ultraviolet behaviour of the solutions of Eqs. (3.1) and 
(4.46) is precisely that of QCD; viz., 



/ 


chiral u, d s c b 




0.0 0.0037 0.082 0.97 4.1 



, r , 9\ P ^^QCD 

M(P) ' (iinrVW" ' ( 8) 



where m is the renormalisation-group-invariant current-quark mass. 

The mass function is a renormalisation group invariant and can be used 
to define a Euclidean constituent-quark mass; viz., 14 



14 The model produces dressed-quarks, which are confined in the sense that their 
propagator does not exhibit a true Minkowski space mass pole; i.e., there is no 
solution of s + M(s) 2 = 0. 
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Fig. 12. Quark mass function obtained by solving the complex formed by Eqs. (3.2), 
(4.44), (4.46), (4.47). NB. A nonzero solution is obtained in the chiral limit. This 
solution is power law suppressed in the ultraviolet [47, 48] and essentially nonper- 
turbative. Thus when the u, d-quark solution melds with the chiral limit solution 
one has entered the domain on which all effects are primarily nonperturbative. The 
transition takes place for m 2 < p 2 < &m 2 g . (Adapted from Ref. [70].) 



{M E f := S \ S = M( S ) 2 . (4.49) 
For the solutions depicted in Fig. 12, one finds (in GeV) 



/ 



Mf 



chiral u,d s c b 

■ (4.50) 

0.42 0.42 0.56 1.57 4.68 v ; 



The ratio Lf := /m/(£) is a measure of the impact of the DCSB mecha- 
nism on a particular flavour of quark. A comparison between Eqs. (4.47) and 
(4.50) shows that for quark flavours with rhf <C m g the effect of DCSB on their 
propagation characteristics is very great V s < &m 2 g . The domain on which the 
impact is large diminishes rapidly as the current-quark mass increases. 

Another way of looking at this is via a constituent-quark cr-term [94, 96] : 

dMf 

°f-=m fiOg^y (4-51) 

It is a renormalisation-group-invariant that can be determined from solutions 
of the gap equation, Eq. (4.46), and is a keen probe of the impact of explicit 
chiral symmetry breaking on the mass function. The ratio af/Mj? plotted 
in Fig. 13 therefore measures the effect of explicit chiral symmetry breaking 
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Fig. 13. Ratio a Q /M§ [Eqs. (4.49) & (4.51)] as a function of current-quark mass. 
It is a measure of the current-quark-mass-dependence of dynamical chiral symme- 
try breaking: a zero value indicates complete dominance of dynamical over explicit 
chiral symmetry breaking; and a value of one, the opposite. The vertical dotted 
lines correspond to the u — d, s, c and b current-quark masses listed in Eq. (4.47). 
(Adapted from Ref. [94].) 

on the dressed-quark mass-function compared with the sum of the effects of 
explicit and dynamical chiral symmetry breaking. It is evident in the fig- 
ure that this ratio vanishes for light-quarks because the magnitude of their 
constituent-mass owes primarily to dynamical chiral symmetry breaking, while 
for heavy-quarks it approaches one. 

An illustration of the chiral limit behaviour can easily be obtained through 
the model in Sect. 2.2, which yields algebraic results in the neighbourhood of 
the chiral limit [97]: M® = ^ G, a = § m, and hence = With 

the parameter values in Ref. [30], this gives an estimate -ff^ f=a 0.04 for u = d 

M o 

quarks. 

For the purpose of a more realistic illustration one may also consider the 
rainbow-truncation defined with Eq. (4.44) and unequal light-quark masses. 
For example, setting m u ((i) = 4MeV and rrid{(i) = 7MeV through Ci = 
1 GeV in Eq. (4.48), one obtains = 0.411 GeV and Mf = 0.416 GeV, and 
cr„ = 0.015 and ad = 0.025. It is interesting to note that the DCSB mechanism 
amplifies the current-quark mass difference in the separation between light- 
quark constituent-quark masses: 

Mf - Mf = 1.5 [m d (Ci) - m„(Ci)] ■ (4.52) 

A value for the vacuum quark condensate, Eq. (4.12), may be obtained 
from the chiral limit solution depicted in Fig. 12 
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-(qq)° c = lim Z 4 ((,A)N c tT D / S, h=0 {q) = (0.275 GeV) 3 . (4.53) 

Since QCD is multiplicatively renormalisable, then 

(?g)?,=Z ra (C,CX?9)g, (4-54) 

where the mass renormalisation constant Z m := Z4/Z2. Applying one- loop 
evolution, Eq. (4.48), to define the vacuum condensate at a typical hadronic 
scale, one therefore obtains [86] 

-(qq}° Ci = (0.24 GeV) 3 . (4.55) 

This condensate might be interpreted as measuring the density of quark- 
antiquark pairs in the vacuum of chiral QCD, in which connection the result 
in Eq. (4.55) corresponds to p qq = 1.8 fm -3 . Were one to assume that this 
vacuum could be viewed as a medium of close-packed spheres, then each would 
occupy a volume V qq = 0.55 fm 3 , which corresponds to a radius r qq = 0.51 fm. 
For comparison, the measured electromagnetic charge radius of a pion can be 
written rv = 0.77r qq and that of a proton, r p = 0.58 r qq . It is therefore clear 
that a veracious understanding of the length-scale defined by the chiral limit 
vacuum quark condensate is a keystone for the bridge between theory and 
experiment. 

The renormalisation-group-improved rainbow-ladder truncation obtained 
with Eq. (4.44) has been employed very successfully. For example, it pre- 
dicted [79, 89] the electromagnetic pion form factor measured [98] at JLab, 15 
and over an illustrative basket of thirty-one calculated quantities, tabulated 
in Ref. [4], its agreement with data has an average relative error of 2% and 
standard-deviation of 15%. 



4.4 Critical mass, DCSB and QCD 

At this point one can return to the question posed at the end of Sect. 2.3; 
namely, whether a critical mass is also associated with DCSB in QCD. This has 
been explored [40] with the renormalisation-group-improved (RGI) rainbow 
truncation of the gap equation's kernel described by Eqs. (4.44), (4.46). 

Equations (3.1) and (4.46) obviously admit the M(p 2 ) = solution in the 
chiral limit, Eq. (3.4). This is the perturbative solution and is analogous to 
Mw in Eq. (2.38). It is moreover plain in Fig. 12 that the chiral limit rainbow 
gap equation also yields a DCSB solution. This capacity is the basis for much 
of the phcnomenological success of the RGI rainbow-ladder truncation. The 
truncation preserves the feature that if M + {p 2 ) = M(p 2 ) > 0, Vp 2 > 0, is 
a solution of the chiral-limit gap equation, then so is M_(p 2 ) := [— M(p 2 )]. 
These solutions are the analogues of M± in Eq. (2.38). 

15 The result is perhaps too good given that pion-loop contributions were omitted. 
Such effects, which appear beyond rainbow-ladder truncation, are important at 
small momentum transfer. For example, they add < 15% to the rainbow-ladder 
result for r„ [99]. 
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Fig. 14. Evolution with current-quark mass, 771(^19), of A(p 2 = 0,£ 2 ) (dimension- 
less), B(p 2 = 0,C 2 ) (GeV) as calculated with uj = 0.4 GeV in Eq. (4.44): solid curve 
- B+(0); dash-dot curve - A+(0); dashed curve - B~(0); dotted curve - A-(0). The 
light dotted curve simply distinguishes the ordinate zero. 



It is apparent in Fig. 14 that this interaction model, too, exhibits a bounded 
domain of current-quark mass on which M_{p 2 ) and M + {p 2 ) exist simultane- 
ously: 

V(m) = {m I < m < m cr } . (4.56) 
The critical current-quark mass depends on u>: 



oj [GeV] 


0.3 0.4 0.5 


m cr [McV] 


71 63 31 


m cr (Ci 9 ) [MeV] 
me, (Ci) [MeV] 


35 31 15 
60 53 26 



where the third and fourth rows report the critical mass at the renormalisation 
scales described above. NB. In most phenomenological applications 0.3 < oj < 
0.4. 

It is noteworthy that when applying Eq. (3.2) only the M±(p 2 ) solutions 
appear; i.e., in this case, for nonzero current-quark mass an analogue of Mw 
in Eq. (2.38) was not found. 

In considering alternative rainbow interaction models, it was determined 
that m cr assumes similar values in most models that provide a reasonable 
description of the same low-energy observables. An exception is the model 
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of Rcf. [30], in which a solution with B(s = 0) < 0, and A(s) and B(s) 
continuous on s <= [0, oo) exists only in the chiral limit. It is also natural to 
enquire after the effect of dressing the quark-gluon vertex. As noted above, 
symmetry ensures that in the chiral limit the gap equation simultaneously 
admits M_(p 2 ) and M + (p 2 ) solutions in this instance, too. The extent of 
the domain of current-quark mass on which the M-(p 2 ) solution persists will 
probably depend on the structure of the vertex. However, that is also true of 
the parameters required to describe a given set of low-energy observablcs. It is 
thus likely that in physical terms the domain T>(rh) will not be much altered 
by realistic vertex dressing. These issues are currently being explored [100]. 

Equation (4.48) is valid for both M + (p 2 ) and M_(p 2 ). In fact, one can 
make a stronger statement, while both M + (p 2 ), M_(p 2 ) exist: 

M + {p 2 ) p2>> ^ M„(j> 2 ). (4.58) 

This follows from asymptotic freedom, which is a feature of the RGI model and 
QCD. One may argue for this result as follows. On the pcrturbative domain 

A (n 2 t 2 \ ~ 1 M±{p2) ~ M ±fr 2 ) (A Ka\ 

A±iP ' C } ~ 1 ' p 2 + Ml(p 2 ) ~ ~ p~ ■ (4 ' 59) 

Hence, the gap equation becomes a single linear integral equation for M±{p 2 ). 
Within the domain on which the preceding steps are valid, that equation can 
be approximated by a linear second-order ordinary differential equation (d.e.) 
(e.g., Ref. [101]). The d.e. is the same for both M±(p 2 ), as is the ultraviolet 
boundary condition, which is determined by the current-quark mass. Thus 
follows Eq. (4.58). 

It is true in general that the model gap equation may usefully be ap- 
proximated by a second-order d.e., which is nonlinear in M(p ) for infrared 
momenta but linear in the ultraviolet [101, 102, 103, 104]. The ultraviolet 
boundary condition for both M±{p 2 ) solutions is still fixed by the current- 
quark mass. However, they have different infrared boundary conditions. The 
d.e. is solved in terms of two independent solutions, which may be described 
as the irregular and regular solutions. The regular solution is power-law sup- 
pressed cf. the irregular form. The irregular solution of the d.e. is the func- 
tional form accessible in perturbation theory. It dominates in the ultraviolet, 
simultaneously fixing both M± (p 2 ) to a value determined by the current-quark 
mass. Any difference between M + {p 2 ) and M_(p 2 ) appears through the regu- 
lar component of the solution, which is essentially nonperturbative and has a 
magnitude that is primarily determined by the infrared boundary condition. 

From Eq. (4.58) and the subsequent discussion it follows that 

fA_ 

a(m(C)) := lim Z 4 (C 2 , A 2 ) N c tr D / S m( %,C), (4.60) 

A ^°° Jq 

S m( %,0 = \ [sf°(qX) - S? (0 (?,C)] , (4-61) 
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Fig. 15. Evolution with current-quark mass, m(^ig), of the massive-quark conden- 
sate defined in Eq. (4.60), calculated with ui = 0.4 GeV and at the renormalisation 
point C = 19 GeV. 

provides a gauge invariant, current-quark-mass-dependent quark condensate 
in QCD, which is finite, well-defined on a bounded interval, and is truly the 
vacuum quark condensate in the chiral limit. The behaviour of (r(m(£)) ob- 
tained with the RGI model employed herein is depicted in Fig. 15. As we saw 
in connection with Fig. 7, here, too, the essentially dynamical component of 
chiral symmetry breaking decreases with increasing current-quark mass. 

It must be stressed that the straightforward definition of a massive-quark 
condensate via the trace of a m ^ dressed-quark propagator is not useful 
because it gives a quantity that is quadratically divergent and therefore very 
difficult to define unambiguously [24, 105]. 

It is now clear that on a bounded interval of current-quark mass, T>(m) = 
{rh | < rh < m cr } with Eq. (4.57), realistic models of QCD's gap equation 
can simultaneously admit two incquivalent dynamical chiral symmetry break- 
ing solutions for the dressed-quark mass- function, M±(p 2 ). These solutions 
arc distinguished by their value at the origin: M + (p 2 ) > and M_(p 2 ) < 0. 
In the ultraviolet they both coincide with the running current-quark mass. 

The pointwise values of both solutions evolve continuously with current- 
quark mass within T>{rh). Beyond the upper boundary, however, only the posi- 
tive solution exists. It is fundamental that M-(p 2 ) is an essentially nonpertur- 
bative quantity whose modification by a current-quark mass can be evaluated 
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Fig. 16. Top line - Integral equation for a planar dressed-quark-gluon vertex. 
"Springs" between dressed-quark lines indicate the dressed-gluon propagator. First 
diagram depicts the rainbow-truncation; second adds a single gluon rung; etc. While 
this series apparently neglects the three-gluon vertex, its effect has satisfactorily been 
modelled via the Ansatz g 2 — > — Cg 2 [61]. The change of sign is an important con- 
sequence of the non-Abelian nature of QCD. Lower two lines - Integral equation for 
a symmetry preserving Bethe-Salpeter kernel that is consistent with the vertex de- 
picted. Beyond rainbow-truncation, symmetry preservation requires that this kernel 
be nonplanar even though the vertex in the gap equation is planar. (Adapted from 
Ref. [109].) 

perturbatively within T>(rh). However, outside that domain such an analysis 
fails because the current-quark mass is large enough to completely destabilise 
this solution. Thus m cr specifics an upper bound on the domain within which 
a perturbative expansion in the current-quark mass can uniformly be valid. 
In the renormalisation-group-improved rainbow-ladder truncation of QCD's 
gap equation described herein the critical current-quark mass corresponds to 
a pseudoscalar meson mass m - ~ 0.45 GeV (mj_ ~ 0.20 GeV 2 ) [106]. 

A value of similar magnitude was deduced in Refs. [107, 108] as the scale 
below which accuracy may be expected from the approximation of observables 
through a perturbative expansion in pseudoscalar meson mass. This scale also 
marks the boundary below which observables should exhibit that curvature 
as a function of pseudoscalar meson mass which is characteristic of chiral 
effective theories. 
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Table 2. Calculated n- and p-meson masses, in GeV, obtained with G — > — C G in 
Eq. (2.22), C = 0.35, and quoted with G = 0.65 GeV in which case m = 0.016G = 
10 MeV. n is the number of interaction rungs retained in dressing the quark-gluon 
vertex, see Fig. (16), and hence the order of the vertex-consistent Bethe-Salpeter 
kernel. NB. n = corresponds to the rainbow-ladder truncation, in which case 
nip = \p2G for m — 0. 





Air 


Air 1 


Air 2 


Mr°° 


7r, m = 














7T, m = 0.011 


0.147 


0.135 


0.139 


0.138 


p, m = 


0.920 


0.648 


0.782 


0.754 


p, m = 0.011 


0.936 


0.667 


0.798 


0.770 



4.5 Beyond rainbow-ladder truncation 

Rainbow-ladder is the leading order in a nonperturbative, systematic and 
symmetry preserving truncation of the DSEs. This is explained and illustrated 
in Refs. [61, 63, 109]. Figure 16 depicts a natural extension of the vertex Ansatz 
in Eq. (2.7). It modifies the gap equation's kernel. In order to preserve the 
axial- vector Ward-Takahashi identity, Eq. (4.2), the Bethe-Salpeter kernel, K, 
must also be modified. That is systematically accomplished for this vertex 
via the procedure indicated in the lower two lines of Fig. 16, as detailed in 
Refs. [61, 109]. 

The model interaction of Eq. (2.22) provides an excellent tool with which 
to illustrate the procedure. When used for the "spring" in the diagrams of 
Fig. 16, it gives the results in Table 2. The first striking observation is that 
to jr = in the chiral limit for every value of n. Plainly, the truncation proce- 
dure preserves the Ward-Takahashi identity, Eq. (4.2). In addition, all of the 
p-TT mass splitting is present in the chiral limit. This answers the questions 
raised on page 3 and in connection with Table 1: the remarkably large dif- 
ference between the ir and p masses owes to DCSB, which forces to be 
unnaturally small. It is also apparent that m n is not very sensitive to the order 
of the truncation. This is another corollary of DCSB for the Goldstone mode; 
namely, the cancellations which ensure = in the chiral limit arc still 
quite effective for small nonzero current-quark masses. Finally, m p does ex- 
hibit some sensitivity to the order of the truncation. With a value of C = 0.35 
the rainbow-ladder truncation is accurate to 21%; the one- loop correction, to 
13%; and the two-loop result, to 4%. 16 Hence, at the two-loop level one is 
certainly very sensitive to details of the model and is thus in a position to 

16 With C — 0.51 the rainbow-ladder truncation overestimates by 29% the complete 
result obtained with the vertex in Fig. 16; with C = 0.25, by 16%; and with 
C = 0.15, by 10%. The discrepancy obviously vanishes as C — > because this 
limit defines the rainbow-ladder truncation. A comparison with these ratios serves 
as a true gauge against which to judge the importance of modifications of the 
vertex in Ref. [61]. 
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current-quark-mass (GeV) 

Fig. 17. Evolution with current-quark mass of the relative difference between the 
qq meson mass calculated in the rainbow-ladder truncation and the exact value 
obtained with C = 0.51. Solid lines: vector meson trajectories; and dashed-lines; 
pseudoscalar meson trajectories. The dotted vertical lines mark the u — d, s, c and 
b current-quark masses. (Adapted from Ref. [61]). 

make a quantitatively accurate map of the interaction between light-quarks 
at long-range via comparison with experiment. These results explain the level 
of accuracy attained with a rainbow-ladder truncation based on Eq. (4.44). 

Reference [61] provides a discussion of the transition between the light- 
quark and heavy-quark sectors, and the relative strength of the corrections 
to rainbow-ladder truncation as this transition is made. Herein the results 
are illustrated with Fig. 17, from which it is apparent that with increasing 
current-quark mass the contributions from nonplanar diagrams and vertex 
corrections are suppressed. Naturally, they must still be included in precision 
spectroscopic calculations. As usual, for small current-quark masses, owing to 
the effects of DCSB, the pseudoscalar channel is a little different. However, 
the trend in this channel becomes the same as that in the vector channel for 
current-quark masses above ~ 2 m a . A contemporary challenge is to generalise 
the procedure described in this section to mesons composed of constituents 
with different current-quark masses. 

5 Baryon Properties 

Significant progress has been made with the study of mesons. While that 
is good, it does not directly impact on the important challenge of baryons. 
Mesons fall within the class of two-body problems. They are the simplest 
bound states for theory. However, the absence of meson targets poses signifi- 
cant difficulties for the experimental verification of predictions such as those 
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reported above. On the other hand, it is relatively straightforward to con- 
struct a proton target but, as a three-body problem in relativistic quantum 
field theory, here the difficulty is for theory. With this problem the current 
expertise is approximately at the level it was for mesons ten years ago; namely, 
model building and phenomenology, making as much use as possible of the 
results and constraints outlined above. 

A natural primary aim is to provide a true picture of the proton's elec- 
tromagnetic form factors and therefrom a reliable determination of the ratio 
depicted in Fig. 1. On the domain of momentum transfer for which there is ap- 
parently a discrepancy between the experimental results; namely, Q 2 > M 2 , 
where M is the nucleon's mass, a veracious understanding of these and 
other contemporary data require a Poincare covariant description of the nu- 
cleoli. This is apparent in applications of relativistic quantum mechanics; 
e.g., Refs. [110, 111, 112, 113, 114]. A different tack follows the formulation 
of a Poincare covariant Faddeev equation [115, 116]. Its foundation is un- 
derstood through the observation that the same interaction which describes 
colour-singlet mesons also generates quark-quark (diquark) correlations in the 
colour-3 (antitriplct) channel [117]. While diquarks do not survive as asymp- 
totic states [61, 63, 109, 118]; i.e., they do not appear in the strong interaction 
spectrum, the attraction between quarks in this channel grounds a picture of 
baryons in which two quarks are always correlated as a colour-3 diquark pseu- 
doparticle, and binding is effected by the iterated exchange of roles between 
the bystander and diquark-participant quarks. 

A first numerical study of this Faddeev equation for the nucleon was re- 
ported in Rcf. [119], and following that there have been numerous more exten- 
sive analyses; e.g., Refs. [120, 121, 122, 123]. It has become apparent that the 
dominant correlations for ground state octet and decuplct baryons are scalar 
and axial-vector diquarks, primarily because the associated mass-scales are 
smaller than the masses of these baryons [124, 125] and the positive parity 
of the correlations matches that of the baryons. Both scalar and axial-vector 
diquarks provide attraction in the Faddeev equation; e.g., a scalar diquark 
alone provides for a bound octet baryon and including axial-vector correla- 
tions reduces that baryon's mass. 

With the retention of axial-vector diquark correlations a quantitative 
description of baryon properties is attainable. Indeed, the formulation of 
Rcf. [121] employs confined quarks, and confined pointlikc-scalar and -axial- 
vector diquark correlations, to obtain a spectrum of octet and decuplet 
baryons in which the rms-deviation between the calculated mass and experi- 
ment is only 2%. The model also reproduces nucleon form factors over a large 
range of momentum transfer [126], and its descriptive success in that appli- 
cation is typical of such Poincare covariant treatments; e.g., Refs. [127, 128, 
129, 130]. 

However, these successes might themselves indicate a flaw in the appli- 
cation of the Faddeev equation to the nucleon. For example, in the context 
of spectroscopy, studies using the Cloudy Bag Model (CBM) [131, 132, 133] 
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indicate that the dressed-nuclcon's mass receives a negative contribution of 
as much as 300-400 McV from pion self-energy corrections; i.e., [134, 135] 
5M + = —300 to —400 McV. Furthermore, a perturbative study, using the 
Faddccv equation, of the mass shift induced by pointlike-7r exchange between 
the quark and diquark constituents of the nucleon obtains [136] SM + = —150 
to — 300 MeV. Unamclioratcd these mutually consistent results would much 
diminish the value of the 2% spectroscopic accuracy obtained using only quark 
and diquark degrees of freedom. 

In addition to masses, pseudoscalar meson loops make important contribu- 
tions to many other baryon properties; e.g., to charge and magnetic radii, and 
magnetic moments [137, 138]. These effects must not be overlooked because 
the size and qualitative impact of meson contributions provide material con- 
straints on the development of a realistic quark-diquark picture of the nucleon, 
and its interpretation and application. 

5.1 Faddeev Equation 

For quarks in the fundamental representation of colour-S'[/(3): 

3 C <g> 3 C <g> 3 C = (3 C © 6 C ) ® 3 C = l c © 8' c © 8 C © 10 c , (5.1) 

and hence any two quarks in a colour-singlet three-quark bound state must 
constitute a relative colour-antitriplet. This fact enables the derivation of a 
Faddccv equation for the bound state contribution to the three quark scat- 
tering kernel [115] because the same kernel that describes mesons so well [4] is 
also attractive for quark-quark scattering in the colour-3 channel. 

In this truncation of the three-body problem the interactions between two 
selected quarks are added to yield a quark-quark scattering matrix, which is 
then approximated as a sum over all possible diquark pseudoparticle terms 
[119]: Dirac-scalar + -axial- vector +[...]. The Faddeev equation thus obtained 
describes the baryon as a composite of a dresscd-quark and nonpointlikc di- 
quark with an iterated exchange of roles between the bystander and diquark- 
participant quarks. The baryon is consequently represented by a Faddeev 
amplitude: 

f = #i + & 2 + *3 , (5.2) 
where the subscript identifies the bystander quark and, e.g., tfq^ are obtained 
from ^3 by a correlated, cyclic permutation of all the quark labels. 

The Faddeev equation is simplified further by retaining only the light- 
est diquark correlations in the representation of the quark-quark scatter- 
ing matrix. A simple Goldstone-theorem-preserving, rainbow-ladder DSE- 
modcl [124] yields the following diquark pseudoparticle masses (isospin sym- 
metry is assumed): 



(qq)jp 


(ud) + 


(ws) + 


(uu)i+ 


(us) 1+ 


m qq (GcV) 


0.74 


0.88 


0.95 


1.05 




(ss)i+ 


(uu)i- 


(us)i- 


(ss)i- 


m qq (GeV) 


1.13 


1.47 


1.53 


1.64 



(5.3) 
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The mass ordering is characteristic and model-independent (cf. Rcfs. [125, 
139], latticc-QCD estimates [140] and studies of the spin-flavour dependence 
of parton distributions [141, 142, 143]), and indicates that a study of the N and 
A must retain at least the scalar [ud] correlations and pscudovector (uu) , {ud) 
and (dd) correlations if it is to be accurate. NB. The spin-3/2 A is inaccessible 
unless pscudovector correlations arc retained. 

The simplest realistic representation of the Faddeev amplitude for the spin- 
and isospin-1/2 nucleon is therefore 

V 3 (p l ,a. l ,n)=Nf +Nf; (5.4) 

namely, a sum of scalar and axial- vector diquark correlations, with (j>i, on, Tj) 
the momentum, spin and isospin labels of the quarks constituting the bound 
state, and P = p\ + P2 +P3 the system's total momentum. 17 For the A, since 
it is not possible to combine an isospin-0 diquark with an isospin-1/2 quark to 
obtain isospin-3/2, the spin- and isospin-3/2 A contains only an axial- vector 
diquark component 

*£(p i ,a un ) = bl + . (5.5) 
The scalar diquark piece in Eq. (5.4) is 

Nf (pi, ai ,n) = [r 0+ (ip [12] ; K)Y a \Z (K) [8(1; P)«(P)]? 3 , (5.6) 

where: the spinor satisfies Eq. (A. 7), and it is also a spinor in isospin space 
with (p + = col(l,0) for the proton and ip~ — col(0, 1) for the neutron; K = 
P1+P2 =■ P{i2},P[i2] =Pi-P2, i := (-p{i2}+2p 3 )/3; A 0+ is a pseudoparticle 
propagator for the scalar diquark formed from quarks 1 and 2, and r° is a 
Bcthc-Salpctcr-likc amplitude describing their relative momentum correlation; 
and S, a 4 x 4 Dirac matrix, describes the relative quark-diquark momentum 
correlation. (S, r° and A° are discussed in Sect. 5.1.) The necessary colour 
antisymmetry of tf^ is implicit in r J , with the Levi-Civita tensor, e clC2C3 , 
expressed via the antisymmetric Gell-Mann matrices; viz., defining 

{H 1 = i\\H 2 = -i\ 5 ,H 3 = iX 2 }, (5.7) 

then e clC2C3 = (H C % 1C2 . [Sec Eqs. (5.32), (5.33).] 
The axial-vector component in Eq. (5.4) is 

N 1+ ( Pliail n) = [t l r^^p ll2] ;K)}ZZ All(K) [K{t-P)u{P)]Z , (5.8) 
where the symmetric isospin-triplet matrices are 



NB. Hereafter isospin symmetry of the strong interaction is assumed; i.e., the u- 
and d-quarks are indistinguishable but for their electric charge. This simplifies 
the form of the Faddeev amplitudes. 



(r°+r 3 ), t c 



r\ t 
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V2 



(r u -r d ) ; 



(5.9) 



and the other elements in Eq. (5.8) are straightforward generalisations of those 
in Eq. (5.6). 

The general form of the Faddeev amplitude for the spin- and isospin-3/2 A 
is complicated. However, isospin symmetry means one can focus on the A ++ 
with it's simple flavour structure, because all the charge states are degenerate, 
and consider 

Dj + = [t+<(ip [12] ;iq]^ 2 All(K) [T> vp (£;P)u p (P)<p + ]ll > ( 5 - 10 ) 

where u p (P) is a Rarita-Schwingcr spinor, Eq. (A. 14). 

The general forms of the matrices S(£; P), AJ,(£; P) and D J/p (£; P), which 
describe the momentum space correlation between the quark and diquark in 
the nucleon and the A, respectively, are described in Ref. [121]. The require- 
ment that S(£; P) represent a positive energy nucleon; namely, that it be an 
eigenfunction of A + (P), Eq. (A. 11), entails 



S(£;P) = Sl (£;P)I D 



n ■ 



£-PId) s 2 (£;P) 



(5.11) 



where £ 2 = 1, P 2 -- 
component, one has 



T. Placing the same constraint on the axial- vector 



At (4 P) = Y, i>n(4 p ) 75 K{£- P) , i = +, 0, • 



(5.12) 



where (£^ = £ 
A}, 



PP V , i-> + 7-M) 



A* —i ip ■ 



1-i-Pv, Al = -iP v , Al 



ri-Al, A 6 V = -£ X -Ai 



(5.13) 



Finally, requiring also that D„ p (£;P) be an eigenfunction of A + (P), one ob- 
tains 

T> vp (£- P) = S A (£; P) 8 up + j 5 A*(t; P) £j , (5.14) 

with and A^ given by obvious analogues of Eqs. (5.11) and (5.12), respec- 
tively. 

It is now possible to write the Faddeev equation satisfied by tf^: 



S(k:P)u(P)] 
K%{k;P)u{P)_ - 

The kernel in Eq. (5.15) is 

M(k,£; P) 



(2*)' 



■M(k,£; P) 



S(£;P)u(P) 
Ai(£; P)u(P) 



M 



on 



(Moi)J 



(Mio)* ( M 



(5.15) 



(5.16) 
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with 

M 00 = r 0+ {k q -£ gq /2;£ qq ) S T (£ qq - k q ) P°\e q ~ k qq /2; -k qq ) S(£ q ) A 0+ (£ gq ) , 

(5.17) 

where: l q = I + P/3, k q = k + P/3, i qq =-£+ 2P/3, k qq = -k + 2P/3, P is 
defined in Eq. (A. 13) and the superscript "T" denotes matrix transpose; and 

(M Ql )i = ^rl + {k q -£ qq l2-£ qq ) 

xs T (e qq - k q ) r 0+ (e q - k qq /2- : -k qq ) s(£ q ) A%(£ qq ) , (5.18) 
(M w yi = r°\k q -i qq /2;e qq ) 

xS T (£ qq - k q ) t 4 Pf(£ q - k qq /2; -k qq ) S(£ q ) A 0+ (£ qq ) , (5.19) 
(M 11 );f„=t-'"r p 1+ (fc g -^/2;^) 
x S T (£ qq - k q ) t* rl\£ q - k qq /2; -k qq ) S{£ q ) Alt (£ qq ) . (5.20) 

The A's Faddeev equation is 

D Ap (fc; P) u p {P) = 4 J ^ M^(k, £; P) D^(£; P) u a {P) , (5.21) 

with 

M^ = t+r^(k q ~£ qq /2;£ qq ) 

x S T (£ qq - k q ) t + rl\£ q - k qq /2- -k qq ) S{Q AH(£ qq ). (5.22) 

The Faddeev equation is illustrated in Fig. 18. It is a linear, homogeneous 
matrix equation whose solution yields the Poincare covariant Faddeev ampli- 
tude, which describes the relative quark-diquark isospin-, linear- and angular- 
momentum-distribution within the baryon. An unamputated Faddeev wave 
function is constructed from the amplitude by attaching the quark and di- 
quark propagators; namely, 

A" + (£ qq )S(£;P)u(P) \ 
A^{£ qq )Al{£-P)u(P)) ■ [ ■ > 

It is this quantity that is most closely related to a wave function in quantum 
mechanics. One can express S(^; P)u(P) in the form of Eq. (5.11). If one labels 
the two functions that appear as ,81,2 then, in the nucleon rest frame, si 2 de- 
scribe, respectively, the upper, lower component of the bound-state nucleon's 
spinor. The lower component is an essentially relativistic effect. It is notewor- 
thy that si has angular momentum zero while s± has angular momentum one. 
The lower component is sizeable even in the Bag Model of the nucleon. It 
follows that the quarks possess a material amount of orbital angular momen- 
tum even in that simple model. Orbital angular momentum is not a Poincare 



~ _ / S(£;P)u(P) \ / 

^-{aWpMp))- 3 ^^ 
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Fig. 18. Pictorial representation of the Faddeev equation, Eq. (5.15). A nucleon 
of four-momentum P is constituted from a dressed-quark (single line, momentum 
p q = k q ) and dressed-diquark (double line, momentum pd = k qq ). Binding is effected 
by an iterated exchange of roles between the bystander and diquark-participant 
quarks, which is described by the kernels M in Eqs. (5.17) - (5.20), (5.22). The 
exchange takes place within the shaded region. 

invariant. However, if absent in a particular frame, it will almost inevitably 
appear in another frame related via a Poincare transformation. Nonzero quark 
orbital angular momentum is the necessary outcome of a Poincare covariant 
description. This is why the covariant Faddeev amplitude is a matrix-valued 
function with a rich structure that, in the baryons' rest frame, corresponds to 
a relativistic wave function with s-wave, p-wave and even ci-wave components. 
(Details can be found in Ref. [144], Sect. 2.4.) 



Faddeev equation kernels 

To complete the Faddeev equations, Eqs. (5.15) & (5.21), one must specify 
the drcsscd-quark propagator, the diquark Bcthc-Salpctcr amplitudes and the 
diquark propagators that appear in the kernels. 

Dressed-quark propagator 

This propagator has the general form given in Eq. (2.3) and herein an overview 
of its properties has already been provided. In solving the Faddeev equation 
described above, one is required to repeatedly evaluate some eight-dimensional 
integrals. Thus, while it is straightforward to obtain a numerical solution 
of a truncation of QCD's gap equation, the utility of an algebraic form for 
S(p) is self-evident. An efficacious parametrisation, which exhibits the features 
described above and has been used extensively in studies of hadron properties, 
is expressed via 

a s (x) = 2mF(2{x + fh 2 )) + F(hx) F(b 3 x) [b + b 2 F{ex)] , (5.24) 
5v{x) = [1 - F(2(x + to 2 ))] , (5.25) 

X + TO 

with x = p 2 /\ 2 , m = to/A, F(x) = (1 — c~ x )/x, as(x) = Xas(p 2 ) and 
av(x) = X 2 ay(p 2 )- The mass-scale, A = 0.566 GeV, and parameter values 18 

18 e = 10 -4 in Eq. (5.24) serves only to decouple the large- and intermediate-p 2 
domains. 
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m b bi b 2 63 



(5.26) 



0.00897 0.131 2.90 0.603 0.185 ' 



were fixed in a least-squares fit to light-meson observables [145, 146]. The 
dimensionless u = d current-quark mass in Eq. (5.26) corresponds to 



and the paramctrisation yields a Euclidean constituent-quark mass, Eq. (4.49), 



which agrees semiquantitatively with that produced by a renormalisation- 
group-improved rainbow-ladder truncation based on Eq. (4.44) [see Eq. (4.50)] . 

Diquark Bethe-Salpeter amplitudes 

The rainbow-ladder DSE truncation yields asymptotic diquark states in the 
strong interaction spectrum. Such states are not observed and their appear- 
ance is an artefact of the truncation. Higher order terms in the quark-quark 
scattering kernel, whose analogue in the quark-antiquark channel do not much 
affect the properties of vector and flavour non-singlct pscudoscalar mesons, en- 
sure that QCD's quark-quark scattering matrix does not exhibit singularities 
which correspond to asymptotic diquark states [61, 63, 109, 118]. Nevertheless, 
studies with kernels that do not produce diquark bound states, do support a 
physical interpretation of the masses, m{ qq ) jP , obtained using the rainbow- 
ladder truncation: the quantity l( qq ) jP = \/m^ qq ^ jP may be interpreted as a 
range over which the diquark correlation can persist inside a baryon. These 
observations motivate the Ansatz for the quark-quark scattering matrix that 
is employed in deriving the Faddeev equation: 



[G qq (k,q;K)]%= £ f J \k; -K) A jP {K) r jP (q; K) , (5.29) 



One practical means of specifying the r J in Eq. (5.29) is to employ the 
solutions of a rainbow-ladder quark-quark BSE. Using the properties of the 
Gell-Mann matrices one finds easily that := r J satisfies exactly 

the same equation as the J~ p colour-singlet meson but for a halving of the 
coupling strength [117]. This makes clear that the interaction in the 3 C (qq) 
channel is strong and attractive. 19 Moreover, it follows as a feature of the 
rainbow-ladder truncation that, independent of the specific form of a model's 
interaction, the calculated masses satisfy 



m = 5.1 MeV, 



(5.27) 



Mia = 0.33 GeV, 



(5.28) 



jp=o+,i+,... 



m iqq)jP > m (qq) _ 



(5.30) 



The same analysis shows the interaction to be strong and repulsive in the 6 C 
(colour-sextet) (qq) channel. 
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This is a useful guide for all but scalar diquark correlations because the part- 
nered mesons in that case are pscudoscalars whose ground state masses, as 
has been shown herein, are constrained to be small by Goldstone's theorem 
and which therefore provide a weak lower bound. For the correlations rele- 
vant herein, models typically give masses (in GeV, recall Eq. (5.3) and the 
associated discussion): 

m {ud) 0+ = 0.74 - 0.82 , m {uu)i+ = m {ud)i+ = m(dd) 1+ = 0.95 - 1.02 . (5.31) 

A solution of the BSE equation requires a simultaneous solution of the 
quark-DSE [125]. However, since in the model under discussion the calcu- 
lations are simplified by parametrising S(p), that expedient should also be 
employed with r J : 

r 0+ (k;K) = J^H a Ci l5 ir 2 F(k 2 /iv 2 0+ ), (5.32) 

t*rf (k; K) = ^H a vy^Ct* F(P/col + ) , (5.33) 

with the normalisation, N J , fixed by 20 

K 2 = -7n 2 p 

(5.34) 
where 

n(K, Q) = tvJ r( q - -K) S(q + Q/2) r(q; K) S T (-q + Q/2). (5.35) 

These Ansatze retain only that single Dirac-amplitude which would represent a 
point particle with the given quantum numbers in a local Lagrangian density: 
they are usually the dominant amplitudes in a solution of the rainbow-ladder 
BSE for the lowest mass J p mesons [79, 86, 88] and diquarks [124, 125]. 

Diquark propagators 

Solving for the quark-quark scattering matrix using the rainbow-ladder trun- 
cation yields free particle propagators for A J in Eq. (5.29). Higher order 
contributions remedy that defect, eliminating asymptotic diquark states from 
the spectrum. The attendant modification of A J can be modelled efficiently 
by simple functions that are free-particle-like at spacclikc momenta but pole- 
free on the timelike axis (see the earlier discussion of confinement, page 24); 
namely, 

20 This is the normalisation condition when the two-body scattering kernel is inde- 
pendent of the system's total momentum. Rainbow-ladder is one such truncation. 
For quark-antiquark bound states the condition has the same form except for the 
replacement S T (-q + Q/2) -> S(q - Q/2) in Eq. (5.35). 



2K„ 



dQ, 



-n(K,Q) 
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W 2 K 2 +), 



A 0+ (K) = 


i 

m l+ 


K( R ) = 


f <W 



K^K V 



1 



F{K 2 fu>l + ),, 



(5.36) 
(5.37) 



'n- 



where the two parameters mjp are diquark pseudoparticle masses and lojf 
are widths characterising r J . It is useful to require additionally that 



dK 2 \m 2 JP 



1 



(5.38) 



K 2 =0 



which is a normalisation that accentuates the free-particle-like propagation 
characteristics of the diquarks within the hadron. 

5.2 Nucleon and A Masses 

The Faddeev equations, Eqs. (5.15) & (5.21), are now completely specified. 
The three-body problem is intrinsically complex and thus it may appear to 
have been a complicated process. However, one has merely combined existing 
information about the one- and two-body sectors of QCD. 



Meson loops and baryon masses 

Before reporting the solution of the Faddeev equation it is worthwhile to 
explicate the effect of pseudoscalar meson loops on baryon masses. A thorough 
discussion is provided in Ref. [122]. One can begin with the leading term in a 
pion-nucleon chiral Lagrangian: 21 



N(x) [i £ - M + ^-757 M t • d^ir(x) + . . .j N(x) . 
The rainbow truncation DSE for the nucleon in this theory is 
- 2 f d 4 k 



E{P) = 3 ? - 



A{k 2 ,ml)$ lb G{P-k)}l lb 



J 4M 2 J (2tt) 4 
wherein the nucleon propagator 

G ( P ) = H wm = : G+ (P) + G~{P) 



P-M-E(P) 



M 



wjv(P) 



^+(P) 



1 



P - lu n (P) +ie 



A_(P) 



1 



Po +lun(P) -is 



(5.39) 



(5.40) 



(5.41) 



(5.42) 



21 That part which describes the pseudoscalar field alone has been neglected. NB. 
This subsection, Eqs. (5.39) - (5.58), employs Minkowski metric for ease of com- 
parison with textbook material. 
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with uj 2 n (P) = P 2 + A/ 2 , and where A±(P) = > ±M)/(2M), P = (w(P),P), 
are the Minkowski space positive and negative energy projection operators, 
respectively; and the pion propagator (w 2 (k) = k 2 + m 2 ) 

A{k\ml) = - (5.43) 



1 



1 



2u> 7r Qt) \_ ko — u) n (k) + is fco + w T (k) — is 



(5.44) 
(5.45) 



As written, the integral in Eq. (5.40) is divergent. It must be regularised 
to give it meaning. In this case the Poincare invariant Paul-Villars procedure 
is useful. It may be effected by modifying the 7r-propagator: 



A(k 2 ,mi) -> A n (k 2 ) = A{k\ m \) + £ c t A(k 2 ,X 2 ), (5.46) 

i=l,2 

and then, with 

A 2 -ml 

c i = - T2 TT > c 2 = T2 yT > ( 5 - 47 ) 

A 2 — A x A 2 — A x 

Eq. (5.46) yields 

^(fc 2 ) = Zi(fc 2 ,m 2 ) J] (A 2 -m 2 )Zi(fc 2 ,A 2 ), 

1=1,2 

(5.48) 

in which case the integrals are convergent for any fixed Ai^- Furthermore, for 
77V Ai — * A 2 = A 

A n {k 2 ) = A{k 2 ,ml) A 2 {k 2 /X 2 , 1) (5.49) 

i.e., the Pauli-Villars regularisation is equivalent to employing a monopolc 
form factor at each irNN vertex: g — > g A(k 2 /X 2 , 1), where k is the pion's 
momentum. Since this procedure modifies the pion propagator it may be in- 
terpreted as expressing compositeness of the pion and regularising its off-shell 
contribution (a related effect is identified in Refs. [147, 148]) but that inter- 
pretation is not unique. 

The contribution to the nucleon's mass from a positive-energy nucleon, 
G + (P), in the loop described by Eq. (5.40) is 

X A M+ = 3g2 [ -f]L_L V A 2 (^y-Af) + 2k 2 (^ Ai + lj n ) 
+ " 16Af 2 J (2tt) 3 w n , = ^ 2 ° l w Xi [u Xi +lu n - M] 

(5.50) 

with u>n = wjy(k 2 ), etc. The connection between this and other mass-shift 
calculations can be made transparent by writing Eq. (5.50) in the form 
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dM +- b7T ^rJ (2^^(k2)K(k2)+ WJV (k2)-A/]' (5 - 51j 

where fjj Nn = g 2 m\ /(16ttM 2 ) and 

k u (k ) := [uj\„ +uj n - M] > ci ■ ■ — , . 

i— 0,1, 2 

(5.52) 

This is useful because for <C Ai — > A2 = A; i.e., on the domain in which 
Eq. (5.49) is valid, one finds algebraically that 

u(k 2 ) = l/(l + k 2 /A 2 ), (5.53) 

which firmly establishes the qualitative equivalence between Eq. (5.50) and 
the calculation, e.g., in Refs. [131, 132, 133, 134]. 

It is instructive to consider Eq. (5.50) further. Suppose that M is very much 
greater than the other scales, then on the domain in which the integrand has 
significant support 

k 2 

Lo N (k 2 )-M^— (5.54) 



2M 



and Eq. (5.50) yields 



i=0,l,2 



It follows that 



d 2 6 A M+ 3g 2 f d 3 k k 2 9 g 2 



(dm 2 ) 2 4M 2 J (2?r) 3 w6(k2) 128tt M 2 

Hence on the domain considered, 



(5.56) 



2 



SaM+ = - J-A_ m ^ + /+ (A l5 A 2 )m 2 + /+(Ai, A 2 ) , (5.57) 

where, as the derivation makes transparent, /(0.1) are scheme-dependent func- 
tions of (only) the regularisation parameters but the first term is regularisation- 
scheme-independent. Given that m 2 oc m, Eq. (4.13), this first term is nonan- 
alytic in the current-quark mass. It is the leading nonanalytic contribution, a 
much touted feature of effective field theory, and its coefficient is fixed by chi- 
ral symmetry and the pattern by which that symmetry is dynamically broken. 
NB. The contribution from G~(P), which produces the so-called Z-diagram, 
is suppressed by 1/M [122]. 

While the leading nonanalytic contribution is model-independent, and thus 
provides a constraint on models that purport to represent QCD, it is not of 
particular quantitative use in this or related studies. The pion and nucleon 
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are both of finite size and hence the regularisation parameters A12, which 
set a compositeness scale for the irNN vertex, must assume soft values; e.g., 
< 600 MeV [126, 128, 149]. Therefore the actual value of the pion-loop contri- 
bution to the nucleon's mass is completely determined by the regularisation- 
scheme-dependent terms. 22 

This is thoroughly explored in Ref. [122], wherein a self-consistent solution 
of the nucleon's gap equation shows that the one-loop result is accurate to 
within 95%. A full consideration leads to the conclusion that the shift in the 
nucleon's mass owing to the nN-loop is (in GeV, for = 1.26): 

-5M+ ~ (0.039 - 0.063) g\ = (0.061 - 0.099) . (5.58) 

Thus one arrives at a robust result: the TrA^-loop reduces the nucleon's mass 
by 100-200 MeV. Extant calculations; e.g., Refs. [131, 132, 133, 134], show 
that the contribution from the analogous 7rZ\-loop is of the same sign and no 
greater in magnitude so that the likely total reduction is 200-400 MeV. Based 
on these same calculations one anticipates that the A mass is also reduced by 
7r loops but by a smaller amount (~ 50- 100 MeV less). 

How is that effect to be incorporated into the quark-diquark picture of 
baryons? Reference [122] argued that it may be included by solving the Fad- 
deev equation with target nucleon and A masses that are inflated to allow 
for the loop corrections. Namely, that a physically sensible picture of the 
quark-diquark piece of the nucleon can be obtained if its parameters are 
chosen not so as to give the experimental masses, but higher values; e.g., 
M N = 0.94 + 0.2 = 1.14 GeV, M A = 1.232 + 0.1 = 1.332 GeV. 

AT and A masses from the Faddeev equation 

The method described in Rcf. [150] is effective for solving Eqs. (5.15) & (5.21). 
Owing to Eq. (5.38), the masses of the scalar and axial-vector diquarks are 
the only variable parameters in the kernels of the Faddeev equations. It is 
natural to choose the axial-vector mass so as to obtain a desired mass for 
the A, and fix the scalar mass subsequently by requiring a particular nucleon 
mass. Two primary parameter sets are presented in Table 3. Set A is obtained 
by requiring a precise fit to the experimental nucleon and A masses, while 
Set B was obtained by fitting to nucleon and A masses that are inflated so 
as to allow for the additional attractive contribution from the pion cloud, as 
described in Sect. 5.2. 

22 From a chiral perturbation theory perspective it might be argued that such A- 
dependent contributions correspond to a short distance effect (d, ~ 1/A ~ 1/3 fm) 
and are therefore absorbed into the a priori unknown counter-term essential 
within that framework. In that approach their effect is thus unobservable. How- 
ever, this is not so in a framework that resolves the internal structure of hadrons. 
In that case the A-dependent contribution from these loops adds unambiguously 
and quantifiably to the quark core defined via the dressed-quark degrees of free- 
dom. The sum predicts a value for the counter-term. 
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Table 3. Mass-scale parameters (in GeV) for the scalar and axial-vector diquark 
correlations, fixed by fitting nucleon and A masses: Set A provides a fit to the 
actual masses; whereas Set B provides masses that are offset to allow for "pion 
cloud" contributions, Sect. 5.2. Also listed is U)jp — rrijp/^/2, which is the width- 
parameter in the (qq)jp Bethe-Salpeter amplitude, Eqs. (5.32) & (5.33): its inverse 
is an indication of the diquark's matter radius. Sets A* and B* illustrate effects 
of omitting the axial- vector diquark correlation: the A cannot be formed and Mjv 
is significantly increased. It is thus plain that the axial- vector diquark provides 
significant attraction in the Faddeev equation's kernel. (Adapted from Ref. [96].) 



set 


M N 


M A 


m + 


m 1 + 


^0+ 




A 


0.94 


1.23 


0.63 


0.84 


0.44=1/(0.45 fm) 


0.59=1/(0.33 fm) 


B 


1.18 


1.33 


0.79 


0.89 


0.56=1/(0.35 fm) 


0.63=1/(0.31 fm) 


A* 


1.15 




0.63 




0.44=1/(0.45 fm) 




B* 


1.46 




0.79 




0.56=1/(0.35 fm) 





It is apparent in Table 3 that a baryon's mass increases with increasing 
diquark mass, and the fitted diquark mass-scales are commensurate with the 
anticipated values, cf. Eq. (5.31), with Set B in better accord. If coupling 
to the axial- vector diquark channel is omitted from Eq. (5.15), then Mjy = 
1.15 GeV and M% ctB = 1.46 GeV, rows labelled A*, B*, respectively. It is thus 
clear that axial- vector diquark correlations provide significant attraction in the 
nucleon. Of course, using the Faddeev equation approach, the A does not exist 
without axial-vector correlations. In Set B the amount of attraction provided 
by axial-vector correlations must be matched by that provided by the pion 
cloud. This highlights the constructive interference between the contribution 
of these two effects to a baryons' mass. It is related and noteworthy that 
m!+ — m + is only a reasonable approximation to Ma — Mn = 0.29 GeV when 
pion cloud effects are ignored: Set A, m 1 + —m + = 0.21 GeV cf. Set B, mi+ — 
m + = 0.10 GeV. Plainly, understanding the N-A mass splitting requires more 
than merely reckoning the mass-scales of constituent degrees of freedom. It 
is curious that for Set B the matter radius, 1/lujp, of the diquarks is smaller 
than for Set A. One might view this as a contraction in the hadron's quark- 
core in response to the presence of a pion cloud that contributes at longer 
range. 

5.3 Nucleon Electromagnetic Form Factors 
Nucleon-photon vertex 

The nucleon's electromagnetic current is given in Eq. (1.6), which introduces 
the Dirac and Pauli form factors, and also the electric and magnetic form 
factors, which are related to the electric-charge-density distribution and the 
magnctic-current-density distribution. In Eq. (1.6), is the nucleon-photon 
vertex. It may be constructed following the systematic procedure of Ref. [151]. 
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Fig. 19. Nucleon-photon vertex which ensures a conserved current for on-shell nu- 
cleons described by the Faddeev amplitudes, S'i,/, calculated in Sect. 5.2. The single 
line represents S(p), the dressed-quark propagator, the double line, the diquark prop- 
agator, and r is the diquark Bethe-Salpeter amplitude, all of which are described 
in Sect. 5.1. Aspects of the remaining vertices are described in Sect. 5.3: the top- 
left image is diagram 1; the top-right, diagram 2; and so on, with the bottom-right 
image, diagram 6. (Adapted from Ref. [152].) 

That approach has the merit of automatically providing a conserved cur- 
rent for on-shell nucleons described by the Faddeev amplitudes which are 
obtained simultaneously with the mass. Moreover, the canonical normalisa- 
tion condition for the nucleons' Faddeev amplitude is equivalent to requiring 
Fi(Q 2 = 0) = 1 for the proton. The vertex has six terms, which are depicted 
in Fig. 19. Hereafter they are briefly described. A full explanation is provided 
in Rcf. [152]. 

Diagram 1 

This represents the photon coupling directly to the bystander quark. It is a 
necessary condition for current conservation that the dresscd-quark-photon 
vertex satisfy the Ward-Takahashi identity: 

Q^ir^x,^) = S-\£ x ) - S-\h) , (5.59) 

where Q = l\ — £2 is the photon momentum flowing into the vertex. Since the 
quark is dressed, Sect. 5.1, the vertex cannot be bare; i.e., ^(^1,^2) In- It 
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can be obtained by solving Eq. (4.25), which was the procedure employed in 
Sec. 4.2. However, since S(p) is parametrised, Ref. [153] can be followed and 
the vertex written [154] 

*T M (*i, * 2 ) = iS A {i\, l\) 7 M + 2fc M [vy ' k A A {1{,1 2 2 ) + A B {£l ; (5.60) 
with k = (£i + £ 2 )/2 and 

S F {llA) = \[F{ll) + F{ll)], A F {£lA) = F[l \Z F i f 2 \ (5-61) 

where F — A,B; viz., the scalar functions in Eq. (2.3). It is critical that in 
Eq. (5.60) satisfies Eq. (5.59) and very useful that it is completely determined 
by the dressed-quark propagator. Moreover, the Ansatz possesses the essential 
feature required by asymptotic freedom that at deep spacelike momenta the 
dressed-quark-photon-vertex evolves to the bare vertex. 

Diagram 2 

This represents the photon coupling directly to a diquark correlation. In the 
case of a scalar correlation, the general form of the diquark-photon vertex is 

rl + (h,£ 2 ) = 2 k, f + (k 2 , k ■ Q, Q 2 ) + /_(fc 2 , k ■ Q, Q 2 ) , (5.62) 
and it must satisfy the Ward-Takahashi identity: 

Q,r° + (£ 1 j 2 ) = n 0+ (e 2 )-n 0+ (t 2 ), n jP (f) - {A jP {i 2 )y\ (5.63) 

The adaption of Eq. (5.60) to this case is 

r° + (£ 1 ,£ 2 ) = k ll A n0+ (£ 2 ,£ 2 ), (5.64) 

with the definition of A n0 + (l 2 ,£%) apparent from Eq. (5.61). Equation (5.64) 
is the minimal Ansatz that: satisfies Eq. (5.63); is completely determined by 
quantities introduced already; and is free of kinematic singularities. It also 
guarantees a valid normalisation of electric charge; viz., 

Inn rl\t',t) = 2£,^n 0+ (£ 2 ) 21, , (5.65) 

owing to Eq. (5.38). NB. The fractional diquark charge has been factored. It 
therefore appears subsequently as a simple multiplicative factor. 

For the case in which the struck diquark correlation is axial-vector, the 
vertex structure is more involved. Nonetheless, there are many constraints 
that may be employed to build a realistic Ansatz. That composed in Ref. [152] 
has two parameters: the magnetic dipole and electric quadrupole moments of 
the axial- vector diquark, and \i+ > respectively. 
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Diagram 3 

This image depicts a photon coupling to the quark that is exchanged as one 
diquark breaks up and another is formed. While this is the first two-loop 
diagram in the current, no new elements appear in its specification: the quark- 
photon vertex was described above. It is noteworthy that the process of quark 
exchange provides the attraction necessary in the Faddeev equation to bind 
the baryon. It also guarantees that the Faddeev amplitude has the correct 
antisymmetry under the exchange of any two dressed-quarks. This key feature 
is absent in models with elementary (noncompositc) diquarks. 

Diagram 4 

This differs from Diagram 2 in expressing the contribution to the nuclcons' 
form factors owing to an electromagnetically induced transition between scalar 
and axial-vector diquarks. The transition vertex is a rank-2 pseudotensor, 
kindred to the matrix element describing the p^*ir° transition [83], and can 
therefore be expressed 

r]%(£i,£ 2 ) = -rj&iliM = ^ r T(e 1 ,£ 2 )e- tapX £ lp e 2 x, (5.66) 

M N 

where 7, a are, respectively, the vector indices of the photon and axial-vector 
diquark. For simplicity, Ref. [152] proceeded under the assumption 

f(£ 1 ,£ 2 ) = Kf ; (5.67) 

viz., a constant, for which a typical value is [126]: 

K t ~2. (5.68) 

In the nucleons' rest frame, a conspicuous piece of the Faddeev wave func- 
tion that describes an axial-vector diquark inside the bound state can be 
characterised as containing a bystander quark whose spin is antiparallel to 
that of the nucleon, with the axial-vector diquark's parallel. The interaction 
pictured in this diagram does not affect the bystander quark but the trans- 
formation of an axial-vector diquark into a scalar effects a flip of the quark 
spin within the correlation. After this transformation, the spin of the nucleon 
must be formed by summing the spin of the bystander quark, which is still 
aligned antiparallel to that of the nucleon, and the orbital angular momen- 
tum between that quark and the scalar diquark. 23 Diagram 4 can therefore 
be expected to impact strongly on the nuclcons' magnetic form factors. 

23 A less prominent component of the amplitude has the bystander quark's spin 
parallel to that of the nucleon while the axial-vector diquark's is antiparallel: this 
© (qq)[ + system has one unit of angular momentum. That momentum is absent 
in the © (qq)o+ system. Other combinations also contribute via Diagram 3 but 
all mediated processes inevitably require a modification of spin and/or angular 
momentum. 
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Diagram 5 and 6 

These two-loop diagrams are the so-called "seagull" terms, which appear as 
partners to Diagram 3 and arise because binding in the nucleons' Faddeev 
equations is effected by the exchange of nonpointlike diquark correlations [151]. 
The new elements in these diagrams are the couplings of a photon to two 
dressed-quarks as they either separate from (Diagram 5) or combine to form 
(Diagram 6) a diquark correlation. As such they are components of the five 
point Schwingcr function which describes the coupling of a photon to the 
quark-quark scattering kernel. This Schwingcr function could be calculated, 
as is evident from the recent computation of analogous Schwingcr functions 
relevant to meson obscrvablcs [93]. However, such a calculation provides valid 
input only when a uniform truncation of the DSEs has been employed to 
calculate each of the elements described hitherto. In the present context, it 
is appropriate instead to employ the algebraic parametrisation of Rcf. [151], 
which is simple, completely determined by the elements introduced already, 
and guarantees current conservation for on-shcll nucleons. 

Calculated results 

In order to place the calculation of baryon observables on the same footing 
as the study of mesons, the proficiency evident in Ref. [4] will need to be ap- 
plied to every line and vertex that appears in Fig. 19. While that is feasible, 
it remains to be done. In the meantime, one can relate a study whose merits 
include a capacity to: explore the potential of the Faddeev equation trunca- 
tion of the baryon three-body problem; and elucidate the role of additional 
correlations, such as those associated with pseudoscalar mesons. 

It is worthwhile to epitomise the input before presenting results. One el- 
ement is the dressed-quark propagator, Sect. 5.1. The form used [145] both 
anticipated and expresses the features that are now known to be true [46, 51]. 
It carries no free parameters, because its behaviour was fixed in analyses of 
meson observables, and is basic to a description of light- and heavy-quark 
mesons that is accurate to better than 10% [70]. 

The nucleon is supposed at heart to be composed of a dressed-quark and 
nonpointlike diquark with binding effected by an iterated exchange of roles 
between the bystander and diquark-participant quarks. The picture is realised 
via a Poincare covariant Faddeev equation, Sect. 5.1, which incorporates scalar 
and axial- vector diquark correlations. There are two parameters, Sect. 5.2: the 
mass-scales associated with these correlations. They were fixed by fitting to 
specified nucleon and A masses. There are no free parameters at this point. 

With the constituents and the bound states' structure defined, only a spec- 
ification of the nucleons' electromagnetic interaction remained. Its formula- 
tion, sketched in Sect. 5.3, is guided almost exclusively by a requirement that 
the nucleon-photon vertex satisfy a Ward-Takahashi identity. Since the scalar 
diquark's electromagnetic properties are readily resolved, the result, Fig. 19, 
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Table 4. Magnetic moments calculated with the diquark mass-scales in Table 3 
for a range of axial- vector-diquark-photon vertex parameters, centred on the point- 
particle values of fj,^ = 2 & Xi+ — li an d Kf = 2, Eq. (5.68). Columns labelled a 
give the percentage-difference from results obtained with the reference parameters. 
Experimental values are: k v := fj, p — 1 = 1.79 &l fi n = —1.91. The magnetic moments 
are expressed in units of nuclear magnetons defined with the calculated nucleon mass. 
(Adapted from Ref. [152].) 

Set A Set B 





Xi+ 






® Kp 














1 


1 


2 


1.79 


-15.3 


1.70 


-5.1 


2.24 


-21.9 


2.00 


-6.2 


2 


1 


2 


2.06 




1.79 




2.63 




2.13 




3 


1 


2 


2.33 


15.4 


1.88 


5.1 


3.02 


21.9 


2.26 


6.1 


2 





2 


2.06 


0.0 


1.79 


0.0 


2.63 


0.0 


2.13 


0.0 


2 


2 


2 


2.06 


0.0 


1.79 


0.0 


2.63 


0.0 


2.13 


0.0 


2 


1 


1 


1.91 


-8.4 


1.64 


-8.4 


2.45 


-10.1 


1.95 


-8.5 


2 


1 


3 


2.21 


8.4 


1.85 


8.3 


2.82 


10.1 


2.31 


8.5 



depends on three parameters that are all tied to properties of the axial- vector 
diquark correlation: Hi+ & Xi+ > respectively, the axial- vector diquarks' mag- 
netic dipole and electric quadrupolc moments; and Kj., the strength of elec- 
tromagnetic axial-vector <-> scalar diquark transitions. Hence, the study of 
Ref. [152] exhibits and interprets the dependence of the nucleons' form factors 
on these three parameters, and also on the nucleons' intrinsic quark structure 
as expressed in the Poincare covariant Faddeev amplitudes. 

Table 4 lists results for the nucleons' magnetic moments, Gjy-(O), where 
N = n,p [152]. It indicates that the moments are insensitive to the axial- vector 
diquarks' quadrupolc moment but react to the diquarks' magnetic moment, 
increasing quickly in magnitude as fii+ increases. As anticipated in connection 
with Eq. (5.68), the nucleons' moments respond strongly to alterations in the 
strength of the scalar «-> axial-vector transition, increasing rapidly 
increased. Set A, which is fitted to the experimental values of Mjv & M^, 
describes the nucleons' moments quite well: k p is 15% too large; and |/z n |, 
16% too small. On the other hand, Set B, which is fitted to baryon masses 
that are inflated so as to make room for pion cloud effects, overestimates k p 
by 47% and \(i n \ by 18%. 

The nucleons' charge and magnetic radii 



r% := -6^1nG^( S ) 



(r%) 2 := -6-lnG&00 



(5.69) 



were also reported in Ref. [152]. The charge radii, particularly that of the 
neutron, were most sensitive to changes in the axial-vector diquarks' electric 
quadrupolc moment, Xi+- That is not surprising given that Xi+ is t ne only 
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Q 2 [GeV 2 ] Q 2 [GeV 2 ] 



Fig. 20. Response of nucleon form factors to variations in the magnetic moment 
of the axial-vector diquark: = 1,2,3; with Xi+ — 1> K f = 2. The legend in 
the top-left panel applies to all; the dashed-line marked by "*" is a fit to experi- 
mental data [155] and the dashed-line marked by "+" in the lower-left panel is the 
fit to G , ] S {Q 2 ) of Ref. [156]; and the horizontal lines in the right panels mark the 
experimental value of the nucleon's magnetic moment. (Adapted from Ref. [152].) 



model parameter which speaks directly of the axial-vector diquarks' electric 
charge distribution. With the point-particle reference values of = 2 & 
Xi+ = 1, and Kf = 2, Set A underestimates the proton radius by 30% and 
the magnitude of the neutron radius by 43%, while for Set B these differ- 
ences are 32% and 50%, respectively. The magnetic radii are insensitive to the 
axial-vector diquarks' quadrupole moment but react to the diquarks' magnetic 
moment as one would anticipate: increasing in magnitude as /!] + increases. 
Moreover, again consistent with expectation, these radii respond to changes 
in Kf , decreasing as this parameter is increased. With the reference parameter 
values both Sets A & B underestimate by approximately 40%. 

Nucleon electromagnetic form factors associated with the tabulated values 
of static properties arc presented in Fig. 20. The figure confirms and augments 
the information in Table 4. Consider, e.g., the electric form factors. One ob- 
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serves that the differences between results obtained with Set A and Set B 
generally outweigh those delivered by variations in the parameters charac- 
terising the axial-vector diquark's electromagnetic properties. The proton's 
electric form factor, in particular, is largely insensitive to these parameters. 
The nucleons' magnetic form factors exhibit the greatest sensitivity to the 
axial-vector diquarks's electromagnetic properties but in this case, too, the 
differences between Set A and Set B arc more si gnificant. For Q 2 > 4 GeV 2 
there is little sensitivity to the diquarks' electromagnetic parameters in any 
curve because the model expresses the diquark current's perturbative limit 
[157]. It is thus apparent from these plots that the behaviour of the nucle- 
ons' form factors is primarily determined by the information encoded in the 
Faddeev amplitudes. 

The results show that the nucleons' electromagnetic properties are sensi- 
tive to the strength of axial-vector diquark correlations in the bound state 
and react to the electromagnetic properties of these correlations. In all cases 
the dependence is readily understood intuitively. However, taken together the 
results indicate that one cannot readily tune the model's parameters to pro- 
vide a uniformly good account of nucleon properties: something more than 
drcssed-quark and -diquark degrees of freedom is required. 



Chiral corrections 



It is appropriate now to examine effects that arise through coupling to pseu- 
doscalar mesons. As with baryon masses, Sect. 5.2, there are two types of con- 
tributions to electromagnetic form factors from meson loops: regularisation- 
schemc-dcpcndent terms, which are analytic in the neighbourhood of rh = 
0; and nonanalytic scheme-independent terms. The leading-order scheme- 
independent contributions for the nucleon static properties presented are [158] 

/ \l-loop ^9a^n _ n v 
\^ p jNA = T Ajr 2f2 m * ' ( 572 ) 

where g A = 1.26, / w = 0.0924 GcV= 1/(2.13 fin), y, v = [i v - [i n . Clearly, the 
radii diverge in the chiral limit, a much touted aspect of chiral corrections. 
While these scheme-independent terms are immutable, at physical values of 
the pseudoscalar meson masses they do not usually provide the dominant 
contribution to observables: that is provided by the regularisation-parameter- 
dependent terms. (We saw this with the baryon masses in Sect. 5.2. See also 
footnote 22.) 

Since regularisation-parameter-dependent parts of the chiral loops are im- 
portant it is sensible to follow Ref. [159] and estimate the corrections using 
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Table 5. Row 1 - static properties calculated with Set B diquark mass es, Tab le 3, 
and = 2, Xi+ = 1, «r = 2: charge radii in fm, with r n := — *J— Jr%); and 
magnetic moments in nuclear magnetons. Row 2 adds the corrections of Eqs. (5.73)- 
(5.75) with A = 0.3 GeV. ? in row n, is the rms relative-difference between the entries 
in rows n and 3. (Adapted from Ref. [152].) 





r p 


r n 


r v 


' n 


M P 






q-(qq) core 


0.595 


0.169 


0.449 


0.449 


3.63 


2.13 


0.39 


+7r-loop correction 


0.762 


0.506 


0.761 


0.761 


3.05 


1.55 


0.23 


experiment 


0.847 


0.336 


0.836 


0.889 


2.79 


1.91 





modified formulae that incorporate a single parameter which mimics the effect 
of regularisation-dcpcndent contributions from the integrals. Thus Eqs. (5.70) 
- (5.72) are rewritten 

/ 2 \l-loop R _ 1 + 5#A i / m x \ fr -70^ 

("Jna =T^^-M^ T ^) J (5-73) 

(K)%-?° PR = - ^| M^fp) + - arctan(A) , 

(5.74) 

/ ^i-ioo P R ^9a m n 2 A 3 . 

= T 7 2f 2 ^ -arctan( — ) , (5.75) 

wherein A is a regularisation mass-scale. NB. The loop contributions vanish 
when the pion mass is much larger than the regularisation scale, as required: 
very massive states must decouple from low-energy phenomena. 

One may now return to the calculated values of the nucleons' static proper- 
ties. Consider the Set B results obtained with = 2, xi+ = 1, Kf = 2. Set B 
was chosen to give inflated values of the nucleon and A masses in order to 
make room for chiral corrections, and therefore one may consistently apply the 
corrections in Eqs. (5.73) - (5.75) to the static properties. With A = 0.3 GeV 
this yields the second row in Table 5: the regularised chiral corrections reduce 
the rms relative-difference significantly. 

This crude analysis, complementing Sect. 5.2, suggests strongly that a ve- 
racious description of baryons can be obtained using dressed-quark and - 
diquark degrees of freedom augmented by a sensibly regulated pscudoscalar 
meson cloud. The inverse of the regularisation parameter is a length-scale that 
may be viewed as a gauge of the distance from a nucleon's centre-of-mass to 
which the pseudoscalar meson cloud penetrates: 1/A ~ § fm is an intuitively 
reasonable value that indicates the cloud is expelled from the nucleon's core 
but materially affects its properties at distances > r^ ™ 3 ; viz., in the vicinity 
of the nucleon's quark-core surface and farther out. 
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Fig. 21. Form factor ratio: /i p G E (Q 2 )/G p lI (Q 2 ). Calculated results: lower band - 
Set A in Table 3; and upper band - Set B. For both, G P E {Q 2 ) was calculated using 
the point-particle values: = 2 & \i+ — 1> an d Kf = 2; i.e., the reference values 
in Table 4. Changes in the axial-vector diquark parameters used to evaluate G P E 
have little effect on the plotted results. The width of the bands reflects the variation 
in G V M with axial-vector diquark parameters. In both cases, the upper border is 
obtained with = 3, Xi+ = 1 an d K f — 2, while the lower has = 1. The data 
are: circles - Ref. [15]; squares - Ref. [16]; and diamonds - Ref. [18]. (Adapted from 
Ref. [152].) 



Form factor ratios 

One is now in a position to return to our discussion of Fig. 1, and in Fig. 21 plot 
the calculated ratio fi p G P E {Q 2 ) / G P M (Q 2 ) . The behaviour of the experimental 
data at small Q 2 is now readily understood. In the neighbourhood of Q 2 = 0, 

M P Q = ° 1 - f [{r P f ~ K) 2 ] , (5.76) 

and because experimentally r p w the ratio varies by less than 10% on 
< Q 2 < 0.6 GeV 2 , if the form factors are approximately dipolc. In the 
calculation described herein, r p > without chiral corrections. Hence the 
ratio must fall immediately with increasing Q 2 . Incorporating pion loops, one 
obtains the results in Row 2 of Table 5, which have r p « r^. The small Q 2 
(long-range) behaviour of this ratio is thus materially affected by the proton's 
pion cloud. 
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It has been emphasised that true pseudoscalar mesons are not pointlikc 
and therefore pion cloud contributions to form factors diminish in magnitude 
with increasing Q 2 . To further exemplify this, it is notable that in a study 
of the jN — > A transition [160], pion cloud contributions to the Ml form 
factor fall from 50% of the total at Q 2 = to < 10% for Q 2 > 2GcV 2 . 
Hence, the evolution of n P G P E {Q 2 ) / G P M {Q 2 ) on Q 2 > 2 GeV 2 is primarily de- 
termined by the quark core of the proton. This is evident in Fig. 21, which 
illustrates that, on Q 2 G (1,5) GeV 2 , [i v G P E (Q 2 )/G P M (Q 2 ) is sensitive to the 
parameters defining the axial- vector-diquark-photon vertex. The response di- 
minishes with increasing Q 2 because the representation of the diquark current 
expresses the pcrturbative limit. 

The behaviour of ^ P G P E (Q 2 )/G P M (Q 2 ) on Q 2 > 2 GeV 2 is determined ci- 
ther by correlations expressed in the Faddeev amplitude, the electromagnetic 
properties of the constituent degrees of freedom, or both. The issue is decided 
by the fact that the magnitude and trend of the results are not materially 
affected by the axial- vector-diquarks' electromagnetic parameters. This ob- 
servation suggests strongly that the ratio's evolution is due primarily to spin- 
isospin correlations in the nucleon's Faddeev amplitude. It is notable that 
while Set A is ruled out by the data, Set B, which anticipates pion cloud 
effects, is in reasonable agreement with both the trend and magnitude of the 
polarisation transfer data [16, 17, 18]. NB. Neither this nor the Roscnbluth [15] 
data played any role in developing the Faddeev equation and nucleon current. 
The agreement between calculation and experiment therefore yields clear un- 
derstanding. (See footnote 3.) 

The calculation is extended to larger Q 2 in Ref. [161], with a prediction 
that the ratio will pass through zero at Q 2 f=a 6.5 GeV 2 . An experiment is 
planned at JLab that will acquire data on this ratio to Q 2 = 9.0 GeV 2 [162]. 
It is expected to begin running around the beginning of 2008. If one adheres 
to a simple interpretation of G P E {Q 2 ) as a Fourier-transform of the electric 
charge distribution within a proton, then the meaning of a zero in this form 
factor is depicted in Fig. 22: it corresponds to a depletion of electric charge 
at the heart of the nucleon. The distribution of magnetic current exhibits 
no such effect. While this picture is not truly valid because the zero appears 
far into the relativistic domain, the parallel it draws between the effects of 
correlations in the wave functions of nuclei and those within the nucleon's 
Faddeev amplitude are useful. 

Figure 23 depicts a weighted ratio of Pauli and Dirac form factors. A 
perturbative QCD analysis [164] that considers effects arising from both the 
proton's leading- and subleading-twist light-cone wave functions, the latter of 
which represents quarks with one unit of orbital angular momentum, suggests 

Q 2 F 2 (Q 2 ) _ constant) q2 >a2) (5 77) 



[lnQ 2 /A 2 ] 2 Fx(Q 2 ) 

where A is a mass-scale that corresponds to an upper-bound on the domain 
of soft momenta. An argument may be made that a judicious estimate of 
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Fig. 22. Three dimensional Fourier transform of a dipole charge distribution com- 
pared with that of a distribution which exhibits a zero; viz., (1 — aQ 2 /m 2 D )/(l + 
Q 2 /m|)) 3 , where ran is the dipole mass and a = 0.1. The zero indicates a deple- 
tion of charge at r = 0, and its relocation to larger r. This profile is familiar from 
nuclear physics: correlations in the proton wave functions of nuclei effect a similar 
redistribution of charge [163]. 

the least-upper-bound on this domain is [152] A = M . The figure hints that 
Eq. (5.77) may be valid for Q 2 > 6GeV 2 . NB. The model for the nucleon 
and its current reviewed herein is consistent with quark counting rules, albeit 
neglecting the anomalous dimensions that arise via renormalisation. However, 
they were also omitted in deriving Eq. (5.77). 

It is worth reiterating that orbital angular momentum is not a Poincare 
invariant. However, if absent in a particular frame, it will in general appear 
in another frame related via a Poincare transformation. Nonzero quark or- 
bital angular momentum is a necessary outcome of a Poincare covariant de- 
scription, and a nuclcon's covariant Faddeev amplitude possesses structures 
that correspond in the rest frame to s-wave, p-wave and even d-wave com- 
ponents. The result in Fig. 23 is not significantly influenced by details of the 
diquarks' electromagnetic properties. Instead, the behaviour is primarily gov- 
erned by correlations expressed in the proton's Faddeev amplitude and, in 
particular, by the amount of intrinsic quark orbital angular momentum [153]. 
NB. This phenomenon is analogous to that observed in connection with the 
pion's electromagnetic form factor. In that instance the so-called axial- vector 
components of the pion's Bethc-Salpcter amplitude, Eqs. (4.9) & (4.10), arc 
responsible for the large Q 2 behaviour of the form factor: they alone ensure 
Q 2 Ft V (Q 2 ) r* constant for truly ultraviolet momenta [65, 66], which is the re- 
sult anticipated from perturbative QCD [165, 166]. These components arc 
required by covariance [64] and signal the presence of quark orbital angular 
momentum in the pseudoscalar pion. 
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Fig. 23. Weighted proton Pauli/Dirac form factor ratio, calculated with A — 
0.94 GeV. The data are described in Fig. 21. The band was calculated using the 
point-particle values: n 1+ = 2 and \i+ = lj an d Kf = 2- Here the upper border is 
obtained with = 1, \i+ = 1 an d K f — 2, and the lower with fj, ± + — 3. (Adapted 
from Ref. [161].) 



5.4 Nucleon weak and strong form factors 

This framework can naturally be applied to calculate weak and strong form 
factors of the nucleon. Preliminary studies of this type are reported in 
Refs. [126, 128]. Such form factors are sensitive to different aspects of quark- 
nuclear physics and should prove useful, e.g., in constraining coupled-channel 
models for medium-energy production reactions on the nucleon. This is im- 
portant to the search for the so-called missing nucleon resonances and the 
related problem of identifying exotic baryons. 25 

For the purpose of illustration of the methods and questions it is worth 
describing briefly the results currently available for three such form factors: 
the axial-vector and pscudoscalar nucleon form factors, which appear in the 

24 This refers to the fact that numerous excited states of the nucleon (N*), which 
are predicted by the SU(6) C2> 0(3) constituent quark model (CQMqxz), are not 
seen in the baryon spectrum that is obtained from amplitude analyses of ttN 
elastic scattering. At this time more than half of the low-lying states predicted 
by CQMqxs are missing [14]. 

25 The "pentaquark" , with a putative valence-quark structure uudds, would be an 
exotic baryon. Confirmation of its existence would likely be a great fillip for 
hadron physics. A perspective on the current status of pentaquark searches is 
provided in Ref. [167]. 
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Fig. 24. Filled circles: <?a(Q 2 ) in Eq. (5.78) calculated in the chiral limit using 
the nucleon Faddeev amplitudes and the axial- vector-nucleon vertex obtained from 
Eqs. (5.82), (5.85) & (5.87). Solid line: dipole fit to the calculation, with mass- 
scale m^, = 1.69 GeV. The shaded band delimits the result's variation subject to 
10% changes in the parameter values in Eq. (5.88). The experimental value of the 
nucleon's axial coupling (<ja ~ 1.27) is marked by a dashed line. (Adapted from 
Ref. [168].) 



axial- vector-nucleon current 

4 M (P',P) = iu(P')^-A 5ll (q;P)u(P) 

= Hp'Hy bn9A(q 2 ) + q^gpiq 2 )] u(P) , (5.78) 

where q = P' — P, j = 1,2, 3 is the isospin index, and the nucleon spinor, 
u(P), is defined in Eq. (A. 7); and the pion-nucleon coupling 

4(P', P) = u(P')Ai(q; P)u(P) = g wNN (q 2 )u(P')ij 5 T j u(P) . (5.79) 

In the chiral limit the pscudovector vertex of Eq. (5.78) takes the following 
form in the neighbourhood of q 2 = [64] 

p ) q = ° rc 8' ular + ^f«AU<r> p ) , (5-80) 

where A^{q; P) is the pion-nucleon vertex and "regular" denotes non-pole 
terms. In addition, g M J| (P',P) = 0. From these observations ensues the 
Goldberger-Treiman relation: 

Mg A (q 2 = 0) = /„ g^ NN (q 2 - 0) , (5.81) 
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Fig. 25. Calculated chiral-limit result for Pa(<2 2 )/(?a(0), solid line, compared with 
data obtained via pion electroproduction in the threshold region, as described in 
Ref. [169]. The band delimits the result's variation subject to 10% changes in the 
parameter values in Eq. (5.88). Dashed line: dipole fit to data with mass-scale = 
1.1 GeV. (Adapted from Ref. [168].) 



where M is the calculated nucleon mass and gA{q 2 ) is solely associated with 
the regular part of the axial-vector vertex. 

The calculation of electromagnetic form factors outlined in Sect. 5.3 sets 
a pattern for determining <7a(<7 2 ), 3p(? 2 ) and g-nNN(q 2 )- To follow this one 
needs to know how a dressed-quark couples to an axial-vector probe. In the 
chiral limit the dressed-quark-axial-vector vertex satisfies Eq. (4.2) with the 
mass-dependent term omitted. Subject to the assumption of isospin symmetry, 
S u = Sd, in which case the chiral-limit axial- vector Ward-Takahashi identity 
is solved by, recall Eq. (5.61), 



rUk;Q) = 75 - ^S A (k 2 + ,k 2 _) 
+ 2k ^ ■ kA A {k\, k 2 _) + 2^ %£ B (fc 2 , fc 2 



(5.82) 



Naturally, Eq. (5.82) is not a unique Ansatz for the dressed-quark-axial-vector 
vertex but it is an adequate starting point. 
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Fig. 26. Filled circles: Chiral limit result for Q 2 gp(Q 2 ) in Eq. (5.78) calculated 
as described in the caption of Fig. 24. Solid line: dipole fit to the calculation, with 
mass-scale mfj = 1.77 GeV. The shaded band delimits the result's variation subject 
to 10% changes in the parameter values in Eq. (5.88). (Adapted from Ref. [168].) 



For the pion-nucleon coupling, one needs the pion's Bcthe-Salpetcr ampli- 
tude and a definition for its extension off pion mass-shell. In chiral QCD one 
has Eqs. (4.8), upon which one can base the Ansatz 

r^(k;Q)=i l5 r^E B (kl,k 2 _), (5.83) 

where Af n is the canonical normalisation constant calculated with this ampli- 
tude via the analogue of Eq. (5.34) (see footnotes 12, 20). 

One also needs to know the following vertices: pion-axial-vector-diquark; 
axial- vector-probe-axial- vector-diquark; and the pion- and axial- vector-probe- 
induccd scalar-diquark «-> axial- vector-diquark transitions. For these one may 
follow Ref. [126]: 

^») 
rf x (jt >p ) 



hi. 



vrl Mr 



(P' ' +P)fiQv ■ 



2M N U 
2 



m M E 
■..7r01 W r\ 

-IK —— Qp . 

Sir 



ff K A01 s tlfl + 2fAr$ ai 



(5.84) 
(5.85) 

(5.86) 
(5.87) 
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Fig. 27. Chiral-limit result for gp(Q 2 ), dash-dot curve. Data obtained via pion 
electroproduction (filled circles) [170] and world average for muon capture at 
Q 2 — 0.88m^ (filled diamond). Dashed curve - current-algebra result; and solid 
curve - next-to-leading order chiral perturbation theory result [169] . (Adapted from 
Ref. [168].) 



where Mq is the Euclidean light-quark constituent-mass, Eq. (5.28), p & p' 
are the incoming and outgoing diquark momenta and Q = (p' — p). Each 
Ansatz introduces one parameter, for which typical values are [126]: 

k" 1 ~ n A1 ~ 4.5 , k* 01 ~ 3.9 , k aq1 ~ 2.1 . (5.88) 

These parameters were used to obtain the results recapitulated herein, with 
the bands representing a variation of ±10%. NB. A scalar diquark does not 
couple to a single pseudoscalar or axial- vector probe. 

With the elements heretofore described one has an analogue for the nu- 
cleoli's chiral current which corresponds to the top four diagrams in Fig. 19. 
This is necessary but not sufficient to guarantee that the axial-vector-nucleon 
vertex automatically fulfills the chiral Ward-Takahashi identity for on-shell 
nuclcons. This situation must be improved and work on that is underway. 

Figure 24 displays the result for the nucleon's axial-vector form factor. 
A comparison with extant data is provided in Fig. 25. The mismatch can 
be attributed to a failure of the axial-vector-nucleon vertex obtained from 
Eqs. (5.82), (5.85) & (5.87) to properly express the diquarks' nonpointlikc 
nature: the momentum evolution thus describes a distribution of axial-charge 
that is too hard. 

In Fig. 26 the result for the nucleon's induced-pseudoscalar form factor is 
depicted. A comparison with data is provided in Fig. 27. The form factor is 



Aspects of Hadron Physics 75 

15i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — I 



10- 



5- 



j i_ 



Q 2 [GeV 2 ] 

Fig. 28. Filled circles: Chiral limit result for (/^]vjv(Q 2 ) in Eq. (5.79) calculated 
using the nucleon's Faddeev amplitude and the irNN vertex constructed from 
Eqs. (5.83), (5.84) & (5.86). Solid line: monopole fit to the calculation, with mass- 
scale mjr = 0.95 GeV. The shaded band delimits the result's variation subject to 
10% changes in the parameter values in Eq. (5.88). The experimental value of the 
nNN coupling (g^NN ~ 13.4) is marked by a dashed line. 



dominated by the pion pole in the neighbourhood of q 2 = —in 2 ,, which for this 
chiral- limit calculation is q 2 ~ 0. The comparison with data is more favourable 
in this instance, particularly once one allows for a shift of the pion pole to 
q 2 = in the chiral-limit calculation. This can be attributed to Eqs. (5.83), 
(5.84) & (5.86); viz., as it is based on Eqs. (4.8), the calculation incorporates a 
fairly accurate representation of pion structure and the pion nucleon coupling. 

This view is supported by the result for g^NN^ 2 ), which is depicted 
in Fig. 28. Within reasonable variation of the parameters that characterise 
the pion-nucleon vertex, the calculated value of g^ NN {0) is consistent with 
standard phenomenology. The result yields a chiral-limit value r® NN ~ 
0.51 ± 0.02 fm. For comparison, a massive-quark value of r^^N ~ 0.3 fm ap- 
pears in Ref. [171], while r^NN ~ 0.93- 1.06 fm is employed in Ref. [172]. 

In order to improve upon these preliminary results, construction must 
be completed of an axial- vector-nucleon vertex that automatically fulfills the 
chiral Ward-Takahashi identity for on-shell nucleons described by the solution 
of the Faddeev equation based on confined-quark and -nonpointlikc-diquark 
degrees of freedom. This will subsequently lead to an improved pion-nucleon 
vertex. In addition, as is known to be necessary for an accurate description of 
nucleon electromagnetic properties, Sect. 5.3, the effect of pseudoscalar meson 
loops on the axial-vector and pseudoscalar couplings must be incorporated. 
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These steps are prerequisites for a veracious description of weak and pionic 
processes. 

6 Epilogue 

Protons and neutrons are the seeds of all the universe's observable matter. 
The standard model of particle physics is supposed to explain their properties. 
However, this theory's perturbative formulation fails spectacularly to account 
for even the simplest bulk properties. Two fundamental, emergent phenomena 
are responsible: confinement and dynamical chiral symmetry breaking. Their 
importance is difficult to overestimate. For example, they determine which 
chemical elements are stable and hence influence even the existence of life. 

Dynamical chiral symmetry breaking (DCSB) is a singularly effective mass 
generating mechanism. It can take the almost massless light-quarks of pertur- 
bative QCD and turn them into the massive constituent-quarks whose mass 
sets the scale which characterises the spectrum of the strong interaction ac- 
cessible at modern experimental facilities. The phenomenon is understood via 
QCD's gap equation, the solution of which delivers a quark mass function 
with a momentum-dependence that connects the perturbative and nonper- 
turbative, constituent-quark domains. 

Despite the fact that light-quarks are made heavy, the mass of the pseu- 
doscalar mesons remains unnaturally small. That, too, owes to DCSB, ex- 
pressed this time in a remarkable relationship between QCD's gap equation 
and those colour singlet Bethe-Salpeter equations which have a pseudoscalar 
projection. Goldstone's theorem is a natural consequence of this connection. 

These features may only be voraciously understood in relativistic quantum 
field theory. They can be viewed as an essential consequence of the presence 
and role of particle-antiparticle pairs in an asymptotically free theory. This 
is apparent via a self-consistent solution of the appropriate Dyson-Schwingcr 
equations (DSEs), which wrap each of QCD's elementary excitations in a 
cloud of virtual particles that is exceedingly dense at low momentum. 

Indeed, the DSEs provide a natural framework for the exploration of 
QCD's emergent phenomena. They are a generating tool for perturbation the- 
ory and thus give a clean connection with processes that are well understood. 
Moreover, they admit a systematic, symmetry preserving and nonperturbative 
truncation scheme, and thereby give access to strong QCD in the continuum. 
On top of this, a quantitative comparison and feedback between DSE and 
lattice-QCD studies is today proving fruitful. 

The existence of a sensible truncation scheme enables the proof of exact 
results using the DSEs. That the truncation scheme is also tractable provides 
a means by which the results may be illustrated, and furthermore a practi- 
cal tool for the prediction of observables that are accessible at contemporary 
experimental facilities. The consequent opportunities for rapid feedback be- 
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tween experiment and theory brings within reach an intuitive understanding 
of nonpcrturbativc strong interaction phenomena. 

Modern, high-luminosity experimental facilities employ large momentum 
transfer reactions to probe the structure of hadrons. They are providing re- 
markable and intriguing new information. For an example one need only look 
so far as the discrepancy between the ratio of electromagnetic proton form fac- 
tors extracted via Roscnbluth separation and that inferred from polarisation 
transfer. This discrepancy is marked for Q 2 > 2 GeV 2 and grows with in- 
creasing Q 2 . At such values of momentum transfer, Q 2 > A/ 2 , where M is the 
nucleon's mass, a true understanding of these and other contemporary data 
require a Poincare covariant description of the nucleon. This can be obtained 
with a Faddeev equation that describes a baryon as composed primarily of a 
quark core, constituted of confined quark and confined diquark correlations, 
but augmented by pseudoscalar meson cloud contributions that are sensed 
by long wavelength probes. Short wavelength probes pierce the cloud, and 
expose spin-isospin correlations and quark orbital angular momentum within 
the baryon. The veracity of the elements in this description makes plain that 
a picture of baryons as a bag of three constituent-quarks in relative s-waves 
is profoundly misleading. 

While there are indications that confinement may be expressed in the 
analyticity properties of the dressed propagators, one cannot say that it is 
understood. Consequently, one pressing task, to which the methods described 
herein can be applied, is the drawing of an accurate map of the confinement 
force between light-quarks within mesons. This will enable a clear connection 
to be established between this force and the realisation of dynamical chiral 
symmetry breaking, and an accounting of the distribution of mass within 
mesons. In addition one may then begin to determine whether the confine- 
ment force-field can be excited to produce exotic systems of light-quarks and 
glue. The same must be done for baryons, but here the difficulties are greater 
because one is confronted with a Poincare covariant three-body problem. New 
and improved tools must be developed, which may extend beyond a rigorous 
grounding of the Faddeev equation. 

It is also crucial to develop the tools necessary for charting the pointwisc 
distribution of quarks and gluons within hadrons. With such in hand, one 
might lay out and apportion the pointwise distribution of mass within the 
hadron, and its evolution with the resolving scale of the probe. That knowl- 
edge could be used to elucidate the impact of confinement and DCSB on 
these distributions; e.g., by exhibiting the effects of quenching these emergent 
phenomena. 

It should now be plain that modern nuclear physics poses numerous chal- 
lenges and may reasonably be expected to provide many new surprises. 
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A Euclidean Space 

It is possible to view the Euclidean formulation of a quantum field theory as 
definitive [54, 173, 174, 175]. That decision is crucial when a consideration 
of nonperturbative effects becomes important. In addition, the discrete lat- 
tice formulation in Euclidean space has allowed some progress to be made in 
attempting to answer existence questions for interacting gauge field theories. 
NB. A lattice formulation is impossible in Minkowski space because the inte- 
grand is not non-negative and hence does not provide a probability measure. 
An heuristic exposition of probability measures in quantum field theory can 
be found in Ref. [176], Chap. 6, while Ref. [54], Chaps. 3 and 6, provides a more 
rigorous discussion in the context of quantum mechanics and quantum field 
theory. 

Our Euclidean conventions are easily made plain. For 4-vectors a, b: 

4 

a • b := a M b v 5^ := 2J a i h , (A.l) 

t=i 

where 8^ v is the Kronecker delta and the metric tensor. Hence, a spacelike 
vector, Qp, has Q 2 > 0. The Dirac matrices are Hermitian and defined by the 
algebra 

{7^,7,1 = 2^. (A.2) 

We use 

7s := - 71727374 , (A. 3) 

so that 

tr [757m7^7p7o-] = -^e^ , £1234 = 1 ■ (A. 4) 

A Dirac-like representation of these matrices is: 





— vt\ 






V IT 


1 


, 74 = 





(A.5) 



where the 2x2 Pauli matrices are: 
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It is possible to derive the Euclidean version of every equation introduced 
above assuming certain analytic properties of the integrands. However, the 
derivations can be sidestepped using the following transcription rules: 

Configuration Space Momentum Space 

1. / d A x M -> -i / d A x E 1. / d A k M I / 



2. ^ -> i 7 B ■ B 2. # -> • fc £ 

3. 4 ^ • A B 3. 4 -> -i^ E ■ A E 

4. -A E ■ B E 4. k u qi* -► -fc £ • ^ 



5. a^a, ->x E -d E 5. fc,,^ -► -fc £ ~ B 



These rules are valid in perturbation theory; i.e., the correct Euclidean space 
integral for a given diagram will be obtained by applying these rules to the 
Minkowski integral. The rules take account of the change of variables and 
Wick rotation of the contour. When one begins with Euclidean space, as we 
do, the reverse is also true. However, for diagrams that represent DSEs which 
involve dressed n-point functions, whose analytic structure is not known a 
priori, the Minkowski space equation obtained using this prescription will have 
the right appearance but it's solutions may bear no relation to the analytic 
continuation of the solution of the Euclidean equation. Any such differences 
will be nonperturbative in origin. It is this fact that makes a choice of metric 
crucial at the outset. 

To return to practical matters, a positive energy spinor satisfies 

u(P, s) (ry • P + M) = = (i 7 • P + M) u(P, s) , (A.7) 

where M is the mass obtained by solving the Faddeev equation and s = ± is 
the spin label. The spinor is normalised: 

u(P, s) u(P, s) = 2M , (A.8) 

and may be expressed explicitly: 




with E = WP 2 +M 2 , 



(A.9) 



( J) • X-=( 1 ) 



For the free-particle spinor, u(P, s) = u(P, 5)^74. 

The spinor can be used to construct a positive energy projection operator: 

A+{P) := 217 S u ( p > s ) s) = 217 H7 ■ P + M) ■ (A ' n) 

s=± 
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A negative energy spinor satisfies 

v(P, s) (i 7 ■ P - M) = = (ij ■ P - M) v(P, s) , (A.12) 

and possesses properties and satisfies constraints obtained via obvious analogy 
with u(P, s). 

A charge-conjugated Bethe-Salpeter amplitude is obtained via 

f(k;P) = C 1 r(-k;P) T C , (A.13) 

where "T" denotes a transposing of all matrix indices and C = 7274 is the 
charge conjugation matrix, = —C. 

In describing the A resonance we employ a Rarita-Schwinger spinor to 
unambiguously represent a covariant spin-3/2 field. The positive energy spinor 
is defined by the following equations: 

(i 7 • P + M) U/i (P; r) = , i^{P- r) = , P^(P; r) = , (A.14) 

where r = —3/2, —1/2, 1/2, 3/2. It is normalised: 

u^P;r') Utl {P;r) - 2M5 r ' r , (A.15) 

and satisfies a completeness relation 

3/2 

— Y, n fi (P;r)u u (P;r)=A + (P)R fiV , (A.16) 
r=-3/2 

where 

I 2 - - 1 - 

= S^Ib - -7m lu + -^PfiPvlD - i-[Pfi7u - Pvln] , ( A -!7) 

with P 2 = — 1, which is very useful in simplifying the positive energy ^\'s 
Faddeev equation. 
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